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Abstract 

We consider the wave and Schrodinger equations on a bounded open connected subset Q, of 
a Riemannian manifold, with Dirichlet, Neumann or Robin boundary conditions whenever its 
boundary is nonempty. We observe the restriction of the solutions to a measurable subset u 
of during a time interval [0, T] with T > 0. It is well known that, if the pair (w, T) satisfies 
the Geometric Control Condition (uj being an open set), then an observability inequality holds 
guaranteeing that the total energy of solutions can be estimated in terms of the energy localized 
in uj x (0,T). 

We address the problem of the optimal location of the observation subset uj among all 
possible subsets of a given measure or volume fraction. We solve it in two different situations. 
First, when a specific choice of the initial data is given and therefore we deal with a particular 
solution, we show that the problem always admits at least one solution that can be regular or 
of fractal type depending on the regularity of the initial data. 

This first problem of finding the optimal uj for each initial datum is a mathematical bench- 
mark but, in view of applications, it is important to define a relevant criterion, not depending 
on the initial conditions and to choose the observation set in an uniform way, independent of 
the data and solutions under consideration. Through spectral decompositions, this leads to 
a second problem which consists of maximizing a spectral functional that can be viewed as 
a measure of eigenfunction concentration. Roughly speaking, the subset uj has to be chosen 
so to maximize the minimal trace of the squares of all eigenfunctions. This spectral crite- 
rion can be obtained and interpreted in two ways: on the one hand, it corresponds to a time 
asymptotic observability constant as the observation time interval tends to infinity, and on 
the other hand, to a randomized version of the deterministic observability inequality. We also 
consider the convexified formulation of the problem. We prove a no-gap result between the 
initial problem and its convexified version, under appropriate quantum ergodicity assumptions 
on fi, and compute the optimal value. 

We also give several examples in which a classical optimal set exists, although, as it hap- 
pens in ID, generically with respect to the manifold SI and the volume fraction, one expects 
relaxation to occur and therefore classical optimal sets not to exist. We then provide spec- 
tral approximations and present some numerical simulations that fully confirm the theoretical 
results in the paper and support our conjectures. 

Our results highlight precise connections between optimal observability issues and quantum 
ergodic properties of the domain under consideration. 
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1 Introduction 



1.1 Presentation of the problems 

Let (M,g) be a smooth n-dimensional Riemannian manifold, n > 1. Let T be a positive real 
number and 51 be an open bounded connected subset of M. In this article we consider both the 
wave equation 

d tt y = A 9 y, (1) 

and the Schrodinger equation 

id t y = A 9 y, (2) 

in (0, T) x ft. Here, A g denotes the usual Laplace-Beltrami operator on M for the metric g. If the 
boundary dfl of f2 is nonempty, then we consider boundary conditions 

By = on (0,T) x <9S1, (3) 

where B can be either: 

• the usual Dirichlet trace operator, By = y\ga, 

• or Neumann, By = f^| Sf2 , where ^ is the outward normal derivative on the boundary dCl, 

• or mixed Dirichlet-Neumann, By = Xr y|an+Xri a^ign' wnere ^ = ToULi with FonTi = 0, 
and XTi is the characteristic function of Ti, i = 0,1, 

• or Robin, By = ^i^ n +(iy\dn^ where fi is a nonnegative bounded measurable function defined 
on dfl, such that J gn (3 > 0. 

Our study encompasses the case where dfl = 0: in this case, ([3]) is unnecessary and f2 is a 
compact connected rt-dimensional Riemannian manifold. The canonical Riemannian volume on M 
is denoted by V gi inducing the canonical measure dV g . Throughout the paper, measurable sets3 
are considered with respect to the measure dV g . 

In the boundaryless or in the Neumann case, the Laplace-Beltrami operator is not invertible 
on Z/ 2 (f2,C) but is invertible in 

Lq(£L, C) = {y € L 2 (Q,, C) | fy(x)dV g = 0}. 

Jn 

In what follows, the notation X stands for the space LQ(n,C) in the boundaryless or in the 
Neumann case and for the space L 2 (£l, C) otherwise. We denote by A = — A g the Laplace operator 
defined on D(A) — {y £ X | Ay £ X and By = 0} with one of the above boundary conditions 
whenever dfl ^ 0. Note that A is a selfadjoint positive operator. For all (y ,^/ 1 ) £ D(A 1 ^ 2 ) x X, 
there exists a unique solution y of the wave equation ([1]) in the space C°(0, T; _D(A 1 / 2 ))nC 1 (0, T; X) 
such that y(0, •) = y°(-) and d t y(0, •) - y 1 ^). 

Let lo be an arbitrary measurable subset of 57 of positive measure. Throughout the paper, the 
notation stands for the characteristic function of lo. The equation ([1]) is said to be observable 
on lj in time T if there exists Cj^\xu) > such that 

4 w) (xc,)ll(y^y 1 )llL(AV 2 ) X x</ T J (l^y(^^)l 2 ^ ^ ^^, (4) 

J If M is the usual Euclidean space K, n then dV g = dx is the usual Lebesgue measure. 
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for all (j/ ,?/ 1 ) £ D{A 1 / 2 ) x X. This is the so-called observability inequality, relevant in inverse 
problems or in control theory because of its dual equivalence with the property of controllability 
(see 48 J. It is well known that within the class of C°° domains f2, this observability property 
holds, roughly, if the pair (w,T) satisfies the so-called Geometric Control Condition (GCC) in Q 
(see 0[I1]), according to which every geodesic ray in Q, and reflected on its boundary according 
to the laws of geometrical optics intersects the observation set uj within time T. In particular, if 
at least one ray does not reach uj within time T then the observability inequality fails because of 
the existence of gaussian beam solutions concentrated along the ray and, therefore, away from the 
observation set. 

A similar observability problem can also be formulated for the Schrodinger equation ([2]) : For 
every y° £ D(A), there exists a unique solution y of |(5J) in the space C°(0, T; D(A)) such that 
y(0, ■) = y°(-). The equation |(5J) is said to be observable on w in time T if there exists (Xui) > 
such that 

c { T s) (x^\\yY D{ A } < / [ \d t y(t,x)\ 2 dV g dt, (5) 

JO Jul 

for every y° £ D(A). It is well known that if there exists T* such that the pair (oj,T*) satisfies 
the Geometric Control Condition then the observability inequality ([5j holds for every T > (see 
[44] ) . Indeed the Schrodinger equation can be viewed as a wave equation with an infinite speed of 
propagation. We refer to [42j for a thorough discussion of the problem of obtaining necessary and 
sufficient conditions ensuring the observability inequality, which is a widely open problem. 

In the sequel, Cj . (Xw) an d C<r(Xw) denote the largest possible nonnegative constants for 
which the inequalities ([4]) and ([5]) hold, that is, 

Cy\x«) = ^{ S J^± X) ^ adt I € D{A 1 / 2 ) x X \ {(0,0)}) , (6) 

and 

- inf ( £ & '*»'';'"' ^ * I / 6 D(A) \ { 0}1 . (7, 

t 112/ Hd(A) J 

They are the so-called observability constants. 

Remark 1. These properties can be formulated in different spaces. For instance, the observability 
inequality ((4]) is equivalent to 

4^(x«)ll(tfV)ll*x(D(^)y < / T / |y(*^)| 2 ^*, (8) 

JO J uj 

for all (y , y 1 ) £ X x (^(A 1 / 2 ))', with the same observability constants. Here the dual is considered 
with respect to the pivot space X. Similarly, the observability inequality (J5j is equivalent to 

4 s) (xu,)||j/°|&< f T l \y(t,x)\ 2 dv g dt, (9) 



for every y° £ X. 

Let (0j)j'eJN* be an orthonormal Hilbertian basis of X consisting of eigenfunctions of A on 17, 
associated with the positiv^l eigenvalues (A 2 ) j£ ]n*. In particular, f n cf>j(x)(f>k(x) dV g is equal to 
whenever j =/= k, and 1 whenever j — k. 

2 Note that, in the Neumann case or in the case dfl = 0, one has X = Lp(n). Otherwise if we would consider 
X = L 2 (f2) in those cases, then we would have Ai = (simple eigenvalue) and ij>i = 1/ iJVg{Q). The fact that in 
those cases we define X = L"q{Q.) permits to keep a uniform presentation for all boundary conditions considered at 
the beginning. 
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Remark 2. Let us provide a spectral characterization of the spaces D(A) and D{A 1 / 2 ). There 
holds 

+00 

D(A) = {y G A | ]T A,% 4>j)h < 
i=i 

and 

+00 

D{A^)={yeX\ ^A|< y ,^)l 2 <+oo}. 
3=1 

In the case of Dirichlet boundary conditions, one hasD(A) = # 2 (fi,C)n.H 1 (fi,C) and DM 1 / 2 ) = 
Hq(Q,C). For Neumann boundary conditions, one has D(A) = {y G _ff 2 (Sl,C) | g^i^Q = 

and / n j/(a;)dFg = 0} and D(A 1 ^ 2 ) = {y G ff^fi.C) | / n j/(ac)dV fl = 0}. In the mixed 
Dirichlet-Neumann case (with T ^ 0), one has D(A) = {y G H 2 (Q,C) | y\r = fl^ = 0} 
and I?(A 1 / 2 ) = ffJ (n,C) = {y G JT^fi.C) | y ]ro = 0} (see e.g. [41]). 

In this article we investigate the two following optimal observability problems. Let L G (0, 1) 
be fixed. 

First problem (optimal design for fixed initial data). 

• Wave equation (JXJ) : given (jAy 1 ) G D(A 1 ^ 2 ) x A, we investigate the problem 
of maximizing the functional 

Gt(x«)= [ T I \d t y(t,x)\ 2 dV g dt, (10) 

JO J UJ 

over all possible measurable subsets oj of VL of measure V g (uS) = LV g (£l), where 
y G C°(0,T;L>(A 1 / 2 ))nC 1 (0,T;X) is the solution of flTJ) such thaty(0,-) = y°(-) 
and §f(0,-) = ?/(•). 

• Schrodinger equation ([2|): given y° G D(A), we investigate the problem of 
maximizing the functional Gt defined by (|10[) over all possible measurable subsets 
uj of SI of measure Vg(u) = LV g (Q), where y G C (0, T; D(A)) is the solution of 
© such that y(0, •) = 

In the analysis of this first problem, the observability inequalities are not required since we 
are dealing with fixed initial data. Accordingly, the optimal set cj, whenever it exists, depends of 
course on the initial data under consideration. As will be shown, this problem is mathematically 
challenging and reveals interesting properties. However it is not relevant enough in view of practical 
applications where the location of the observation or sensors is expected to be uniform with respect 
to the data and solutions under consideration. 

Consequently, we introduce the following second problem, of a spectral nature, in which, to 
some extent, all possible solutions are taken into consideration in the optimality criterion. 

Second problem (uniform optimal design) We investigate the problem of maxi- 
mizing the spectral functional 

J( XU )= >f f ^(xfdVg, (11) 

over all possible subsets cuofQof measure V g (u>) — LV g (Q). 
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A relevant and natural criterion would certainly consist in maximizing the observability constant 
over all possible subsets uj of ft of measure V g {oS) = LV g (fl) for a given time T > 0. Settled as 
such this problem is however very difficult to handle. Indeed, using an Hilbertian expansion of 
the solutions of ([T]) or ([2j in the basis of the eigenfunctions of the Laplacian operator, this leads 
to inequalities in which the presence of crossed terms makes it difficult to analyze the existence 
and possible nature of the optimal sets. Furthermore, this criterion depends on the time interval 
[0, T) while the spectral one above is independent of T and is of diagonal nature, not involving any 
crossed term. 

The difficulty related with the cross terms already appears in one-dimensional problems (see 
[54]). Actually, this question is very much related with classical problems in non harmonic Fourier 
analysis, such as the one of determining the best constants in Ingham's inequalities (see [551 136] 1. 

In Section l2~2l we describe how the spectral criterion defined above can be derived by various 
averaging processes applied to the original problem of optimizing the observability constant. The 
first one is to perform a time averaging process, leading to interpret the criterion J(xu) defined 
by ([TT]) as a time asymptotic observability constant as T tends to +00. The second one consists 
of randomizing the initial data of the wave or Schrddinger equation under consideration, which 
leads to interpret J{xu>) a s a randomized observability constant, corresponding to a randomized 
observability inequality (see Section 12721 for details). 

These notions of time asymptotic or randomized observability inequalities are new and happen 
to be better fitted to provide a relevant answer to the problem of optimal observability. We 
provide in Section T2.2| precise relations between these new observability constants and their classical 
deterministic versions. 

Note that, when the spectrum of A is not simple, this spectral second problem depends a 
priori on the choice of the orthonormal basis (^j)iew °f eigenfunctions. When the spectrum is 
not simple it is natural to consider an intrinsic variant of the second problem by considering the 
infimum over all possible normalized eigenfunctions (see Section 14.81) . 

1.2 Brief state of the art 

The literature on optimal observation or sensor location problems is abundant in engineering appli- 
cations (see e.g. [39, 51, 60]|63ll65] and references therein), but very few mathematical theoretical 
contributions do exist. In engineering applications, the aim is to optimize the number, the place 
and the type of sensors in order to improve the estimation of the state of the system. Fields of 
applications are very numerous and concern for example active structural acoustics, piezoelectric 
actuators, vibration control in mechanical structures, damage detection and chemical reactions, 
just to name a few of them. In most of these applications however the method consists in ap- 
proximating appropriately the problem by selecting a finite number of possible optimal candidates 
and of recasting the problem as a finite dimensional combinatorial optimization problem. Among 
these approaches, the closest one to ours consists of considering truncations of Fourier expansion 
representations. Adopting such a Fourier point of view, the authors of [30l [31] studied optimal 
stabilization issues of the one-dimensional wave equation and, up to our knowledge, these are the 
first articles in which one can find rigorous mathematical arguments and proofs to characterize the 
optimal set whenever it exists, for the problem of determining the best possible shape and position 
of the damping subdomain of a given measure. In [5] the authors investigate the problem mod- 
eled in [60] of finding the best possible distributions of two materials (with different elastic Young 
modulus and different density) in a rod in order to minimize the vibration energy in the structure. 
For this optimal design problem in wave propagation, the authors of [8] prove existence results and 
provide convexification and optimality conditions. The authors of [1] also propose a convexification 
formulation of eigenfrequency optimization problems applied to optimal design. In |24j the authors 
discuss several possible criteria for optimizing the damping of abstract wave equations in Hilbcrt 
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spaces, and derive optimality conditions for a certain criterion related to a Lyapunov equation. In 
[54] we investigated the second problem presented previously in the one-dimensional case. We also 
quote the article |55j where we study the related problem of finding the optimal location of the 
support of the control for the one-dimensional wave equation. 

1.3 Short description of the main results of this article 

In this article we provide a complete mathematical analysis of the two optimal observability prob- 
lems settled in Section [O] The article is structured as follows. 

Section [5] is devoted to spectral considerations and to state and prove results interpreting the 
second problem in terms of a time averaged or a randomized observability inequality (see Corollary 
[T] and Theorem [T]) . In particular it is shown how the time averaging or the randomization with 
respect to initial data permit to rule out crossed terms and lead to the spectral criterion (TTTj) 
considered in the second problem. 

In Section [3j we solve the first problem, that is the optimal design problem for fixed initial 
data. The main result of this section is Theorem [3j which provides a sufficient condition ensuring 
existence and uniqueness of a solution of the first problem (see the more precise statement in 
Section [53]). 

Theorem. // the initial data under consideration belong to a suitable class of analytic functions, 
then the first problem has a unique solution u>, which has a finite number of connected components. 

Here, the optimal set w is unique up to some subset of zero measure. Proposition [T] (Section 
I3.2[) shows that the above sufficient condition is, in some sense, sharp, since there exist initial data 
of class C°° such that the first problem has a unique solution lj, which is a fractal set and thus 
has an infinite number of connected components. An explicit example is built in Appendix |A"1 

In Section |4j we focus on the second problem (uniform optimal design), which is of a spectral 
nature and, thus, independent on the initial data. As proved in Section [5] this problem corresponds 
to computing the maximal possible value of the time asymptotic or of the randomized observability 
constant. We first provide in Section l4.1l a convexified version of the problem, by considering the 
convex closure of the set 

Ul = {Xw | w is a measurable subset of ft of measure V 3 (uj) = LV g (fl)} 

for the L°° weak star topology, that is 

U L = {ae L°°(n,[Q,l]) | f a(x)dV g = LV g (Q)}. 

The convexified second problem then consists of maximizing the functional 

J ( a ) = M, / a{x)4>3{xf dV g 

over U l ■ Our main results are Theorems |4] and [5] (stated in Section 14. 2|) , whose contents are 
roughly the following. 

Theorem. 1. Assume that there exists a subsequence of the sequence of probability measures 
fij — cf)j dV g converging vaguely to the uniform measure v ^ dV g (Weak Quantum Ergodicity 
assumption), and that the sequence of eigenf unctions <f>j is uniformly bounded in L°°(f2). 
Then 

sup inf / d>j(x) 2 dV„ = sup inf / a(x)d)j[x) dV„ = L, 
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for every L g (0,1). In other words, there is no gap between the second problem and its 
convexified version. 

2. Assume that the whole sequence of probability measures fij — 4> 2 dV g converges vaguely to the 
uniform measure y-jj^j dV g (Quantum Unique Ergodicity assumption), and that the sequence 
of eigenfunctions <pj is uniformly bounded in L' 2p (fl), for some p G (1, +oo]. Then the 
supremum of J over the subset of Ul of all characteristic functions of Jordan measurable 
subsets is as well equal to L. 

The assumptions of the above result are sufficient but not necessary to derive such a no-gap 
statement, as it is shown when 51 is a two-dimensional disk with Dirichlet boundary conditions 
(see Proposition^, in spite of the fact that the eigenfunctions do not equidistribute as the eigen- 
frequencies increase, as illustrated by the well known whispering galleries effect (see discussion in 
Section O. 

In Section 14.31 we comment on these quantum ergodicity assumptions, that are well known 
in mathematical physics to be related with concentration phenomena of eigenfunctions (scarring 
phenomena). We briefly survey the main results of this theory and show their intimate relations 
with the optimization problems under consideration in the present article. 

In Section r4.4l we provide some results on the existence of an optimal set achieving the supremum 
in the above problem, and formulate some open problems. 

Note that Theorems [4] and [5] cannot be inferred from usual T-convergence results. Theorem [4] 
is proved in Section l4~5l The proof of Theorem [5] in Section l4~6l is of a completely different nature. 
It is a constructive proof of a maximizing sequence of optimal sets which permits to establish that 
it is possible to increase the values of J by considering subsets having an increasing number of 
connected components. 

In Section 14.81 we define and study an intrinsic version of our second problem, which does 
not depend on the choice of an orthonormal basis of eigenfunctions. We consider the problem of 
maximizing the quantity inf^ e £ J 4>{x) 2 dV g over Ul, where £ denotes the set of all normalized 
eigenfunctions of A. For this problem we have a result similar to the one above (Theorems [6] and 
[7]) , and moreover in this intrinsic problem we are able to provide an explicit example where a gap 
occurs between the problem and its convexified formulation (Proposition [3]), by considering the 
unit sphere in M 3 or the unit half-sphere with Dirichlet boundary conditions, and certain quantum 
limits of a Dirac type. 

These results constitute the main contributions of this article. They show precise connections 
between optimal observability issues and quantum ergodicity properties of f2. Such a relation was 
suggested in the early work [2U] concerning the exponential decay properties of dissipative wave 
equations. 

Section [5] is devoted to the study of a finite-dimensional spectral approximation of the second 
problem, namely the problem of maximizing the functional 

Jn(Xu>) = min / (j)j{xfdV g 

over Ul- In Theorem |8] we derive a T-convergence property of Jn towards J for the weak star 
topology of L°°. In particular, the sets optimizing Jjv constitute a maximizing sequence for the 
convexified version of the maximization problem for J, and this, without geometric or ergodic- 
ity assumptions on Q. Of course, then, under the assumptions of the above theorem, these sets 
constitute a maximizing sequence for the original optimization problem as well, without convexifi- 
cation. We also prove the existence and uniqueness of an optimal set lu n , which is shown to have 
a finite number of connected components for every integer N. We also present several numerical 
simulations showing the shapes and complexity of these sets. 
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In Section [51 we provide further comments. First, in Section 16.11 for the second problem 
(|45|) we analyze possible ways of restricting the classes of domains under consideration to ensure 
compactness properties. But then, of course, the maximal value of J diminishes. In Section [6. 2 1 we 
extend our main results (for the second problem) to a natural variant of observability inequality for 
Neumann boundary conditions or in the boundaryless case. Section [631 is devoted to the analysis of 
a variant of observability inequality for Dirichlet, mixed Dirichlet-Ncumann and Robin boundary 
conditions, involving a H 1 norm. We show that the time averaging or the randomization lead to a 
slightly different spectral criterion (Theorem [TU)) . In contrast to the previous results, in this new 
situation the uniform optimal design problem has an optimal solution whenever the volume fraction 
is large enough, the optimal set being determined by taking only into account a finite number of low 
frequency modes (see Theorem [TT] for a precise statement). Numerical simulations illustrate this 
result. Finally, in Section [63] we show that the problem of maximizing the observability constant, 
studied along the paper, is by a classical duality argument equivalent to the optimal design of 
the control problem of determining the optimal location of internal controllers, for the wave and 
Schrodinger equations. 



2 Preliminaries: spectral considerations 

2.1 Spectral expansion of the functional Gt(Xoj) 

We recall that we have fixed a Hilbertian basis {4'j)je'iN* of X consisting of eigenfunctions of A, 
associated with the real eigenvalues (^j)jeiN* ■ Using a series expansion of the solutions of the wave 
or Schrodinger equation in this Hilbertian basis, our objective is to write the functional Gt defined 
by f| lOj) in a more suitable way for our mathematical analysis. 

Wave equation Q. For all initial data (y°, y 1 ) E D(A 1 / 2 )xX, the solution?/ e C°(0, T; D(A 1 ^ 2 ))n 
C X {Q,T\X) of (ID such that y(0, ■) = y°{] and d t y(0, ■) = y x (') can be expanded as 

-j-oo 

y(t, x)=J2 foe"'* + b^*) <Pj(x), (12) 
j=i 

where the sequences (AjOj^gN* and (Xjbj)j^* belong to £ 2 (<C) and are determined in terms of 
the initial data (y ,?/ 1 ) by 



' / ' y°( x )^(x)dV g --^- I y 1 (x)cf> j (x)dV g 



i 



9 ' 



h = l[ I q y°(x)M*) dV g + -J^ y'ix^ix) dV g 



(13) 



for every j G IN*. Moreover, 



+00 

ua/V)ii V/ 2 )xx = 2 E A KM 2 + n 2 )- ( 14 ) 
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Plugging Unj) into {TDJ leads to 



Jlj 

+00 



+00 



x i ( a J 



AX* t 



i\jt\ 



dV g dt 



T^~> r. 



9' 



where 



Ojfc = / ( a i e 



)(a fc e" 



b k e iX ^)dt. 



(15) 



(16) 



The coefficients a^-fe, (j, A;) € (IN*) 2 , depend only on the initial data (y ,?/ 1 ), and their precise 
expression is given by 



2a,a k . 



Xi — X 
'i 



2bjak 



T 



(Ai-A fc )- e ^-A*)T 



Aj + Afe 
whenever Aj 7^ Afe, and 



(A, + A fc) ^ ) e^+^ + sin ( (A, - A*)± ) 



2a,jbk . 
- — J —— sir 

Aj + Afc 

26,-6. 



(Aj+A fe )- ) e ^-+^)i 



(17) 



a^fe = T(ajdk + bjbk) 



X< - X 



(18) 



when A,- = A* 



Remark 3. In dimension one, consider SI = [0, n] with Dirichlet boundary conditions. Then 
4>j{x) — ■\f^sin(jx) and Xj = j for every j £ IN*. In this one-dimensional case, it can be noticed 
that when the time T is a multiple of 27r all nondiagonal terms vanish. Indeed, if T = 2p?t with 
p £ IN* , then = whenever i ^ j, and 

a iJ =p7r(|o i | 2 + |6 i | 2 ), (19) 

for all € (IN*) 2 , and therefore 



+00 . 
G2 P n(Xu) = ^ A 2 ajj / sin 2 (jx)dx. 



(20) 



Hence in that case the functional G^pw does not involve any crossed terms. The second problem 
for this one-dimensional case was studied in detail in [54] . 

Schrodinger equation ©. For every y° £ D(A), the solution y e C°(0,T;D(A)) of © such 
that 2/(0, •) = 2/°( - ) can be expanded as 



+00 



where the sequence (A 2 Cj) je ]N* belongs to ^ 2 (C) and is determined in terms of y° by 



y°(x)4> j (x)dV g 



(21) 



(22) 
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for every j G IN*. Moreover, 

+00 

lly°ll^) = E A !N 2 - ( 23 ) 

Plugging dH]) into dTUJ) leads to 

2 

„T „ +00 +00 . 

G T (Xu)= / / ^Xfee^faix) dV g dt= A|A|a ife / ^(a:)&(a;)dV^ (24) 



with 



~c k I e^>« dt = |^ sin ((A? - Al)|) e^~^ , (25) 



whenever j ^ k, and = |cj| 2 T whenever j = k. 

2.2 Two motivations for studying the second problem 



In this section, as announced at the end of Section 11.11 we provide two motivations for studying 
the second problem, consisting of maximizing the functional J defined by (fTl"]) over the set Ul . 

First motivation: averaging in time / time asymptotic observability constant. 

First of all, we claim that, for all (y°, y 1 ) G D{A 1 / 2 ) x X, the quantity 



1 " T 



TJo 



\d t y(t,x)\ 2 dV g dt, 



where y G C°(0, T; D(A 1 / 2 )) n C^OjTjX) is the solution of the wave equation Q such that 
y(0, •) = y°{-) and d t y(0, •) = y 1 ^), has a limit as T tends to +00. We refer to lemmas Q] and [2] for 
a proof of this fact. This leads to define the concept of time asymptotic observability constant 

^W = wfh ^ij^t X)?dV9dt I (yV)e^ 1/ V*\{(0,0)}). (26) 

This constant appears as the largest possible nonnegative constant for which the time asymptotic 
observability inequality 

c { J\x.)\\{y\y l )\\l {A v* )xX ^Km I f [ \d t y[t,xf\dv g dt, (27) 

holds for all y°{-) G Z?(A 1 / 2 ) and j/ x (-) G X. 

Similarly, for the Schrodinger equation, we define 

- inf ( T lim I ,0 g UQ} \ (2g) 

This constant is the largest possible nonnegative constant for which the time asymptotic observ- 
ability inequality 

Cg>(Xu)\\v°\\D<.A) < t 1 ™^ C i \W>*?\W B to, (29) 



Jul 



holds for every y°(-) G D(A). 
We have the following results. 
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Theorem 1. For every measurable subset ui of CI, there holds 



x 



(X.) = C^fc.) = inf ' u ~™ - | (c,W e *»(C) \ {0} , 



2 



ESN 

where U is the set of all distinct eigenvalues A& and /(A) = {j G IV* | Aj = A}. 
Corollary 1. There holds 

for every measurable subset ujoffl. 

If the domain 57 is such that every eigenvalue of A is simple, then 

2CL ff) (x.) = C«?(x u ) = inf / ^(x) 2 dV a = J(x„), 

for every measurable subset lo of 57. 

The proof of these results arc done in Section [2751 Note that, as is well known, the assumption 
of the simplicity of the spectrum of the Dirichlet-Laplacian is generic with respect to the domain 57 
(see e.g. [50, 64, 33 ). The spectrum of the Neumann-Laplacian is also known to consist of simple 
eigenvalues for many choices of 57. For instance, it is proved in |33j that this property holds for 
almost every polygon of M 2 having N vertices. 

Remark 4. It follows obviously from the definitions of the observability constants that 

hmsup ^l ^ <C4 W) (XJ and limsup ^ c£>( X »), 

for every measurable subset uj of 57. However, the equalities do not hold in general. Indeed, 
consider a set 57 with a smooth boundary, and a pair (w,T) not satisfying the Geometric Control 

Condition. Then there must hold (x«) = 0- Besides, J{Xu) may be positive. 
An example of such a situation for the wave equation is provided by considering 57 = [0,7r] 2 
with Dirichlet boundary conditions and L = 1/2. It is indeed proved further (see Lemma [6] 
and Remark I22[) that the domain u) = {(x, y) € 57 | x ^ tt/2} maximizes J over Ul, and that 
J(Xu>) = 1/2- Clearly, such a domain does not satisfy the Geometric Control Condition, and one 
has Cf = 0, whereas CgP(xu) = 1/4- 

Another class of examples for the wave equation is provided by the Bunimovich stadium (shaped 
at the top right of Figure [2] further) with Dirichlet boundary conditions. Setting 57 = R U W, 
where R is the rectangular part and W the circular wings, it is proved in [18j [19] that, for any 
open neighborhood lo of the closure of W (or even, any neighborhood lj of the vertical intervals 
between R and W) in 57, there exists c > such that J 4>j( x ) 2 dx ^ c for every j € IN*. It follows 

that J(xu) > 0, whereas . (xo») — since lj does not satisfy the Geometric Control Condition. 
It can be noted that the result still holds if one replaces the wings W by any other manifold glued 
along R, so that 57 is a partially rectangular domain. 

We are not aware of such kinds of examples for the Schrodinger equation, although there exist 
some configurations for which Cj? (Xu) = 0. For instance, for 57 = 5 2 , the unit Euclidean sphere 
of M 3 , it is well known (see for instance [37] and Remark [T9l further) that, if (^^jev is the 
usual orthonormal basis of spherical harmonics, then a subsequence of 2 converges to the Dirac 
measure along the equator. Therefore, if a subset a; of 57 does not contain any neighborhood of this 
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equator, then C^'(Xu) — Co£ ; (x w ) = J(Xu) — (see Corollary [JJ) . Note that the same situation 
occurs in the unit disk of the Euclidean plane, choosing any subset oj compactly included in the 
disk, since there exists a subsequence of the squares of the usual Dirichlet-Laplacian eigenfunctions 
concentrating on the boundary of the disk (see e.g. [40], see also Section 14771 further) . 

Second motivation: averaging with respect to initial data / randomized observability 
constant. 

The observability constants in and ([7]) are defined as an infimum over all possible (deterministic) 
initial data. We are going to modify slightly this definition by randomizing the initial data in 
some precise sense, and considering an averaged version of the observability inequality with a new 
(randomized) observability constant. To make this point precise, we consider spectral expansions 
of the solutions of the wave and Schrodinger equations, and we get 



— inf 

2 (A;;aj),(Aj6 3 )e^(C) ./o 
E+~A 3 2 (|a J | 2 + |6 3 | 2 ) = l 



A\j t 



bje 



-iXj t 



) 4>j{x) 



dV g dt, 



and 



C ( t\xu,)= inf 

(x 2 jCj )ee 2 (C) Jo J u 



dV g dt. 



The coefficients a,j, bj and Cj in the expressions above are the Fourier coefficients of the initial 
data, defined by (fl~3|) and ([22]) respectively. 

Following the works of N. Burq and N. Tzvetkov on nonlinear partial differential equations 
with random initial data (see [T2l [15j [TBI U-2]) using early ideas of Paley and Zygmund (see [52]). 
we randomize these coefficients by multiplying each of them by some well chosen random law. 
This random selection of all possible initial data for the wave equation ([JJ) consists of replacing 
C^ V \xoj) by the randomized version 



a 



(W) 
T,rand 



£+3 A?(K-| 2 +|f>il 2 )=i 



E 



/ 2 \ 

/ / E^^^'-W^^w dV s dt 

(30) 

where t j)jeJSi* and ifi% j)jevi* are two sequences of independent Bernoulli random variables on 
a probability space (X,A,P), satisfying 

Vffij = ±1) = HPh = ±1) = \ and WljPlk) = 

for every j and k in IN* and every v G X. Here, the notation E stands for the expectation over 
the space X with respect to the probability measure P. In other words, instead of considering the 
deterministic observability inequality (U]) for the wave equation ([TJ, we consider the randomized 
observability inequality 



a 



(W) 
7\rand 



(X^)ll(2/°,y 1 )llD( J 4i/2)xX 



s; E 




\d t y v (t,x)\ 2 dV g dt 



(31) 



for all y°(-) G D{A 1 / 2 ) and y x { ) G X, where y v denotes the solution of the wave equation with 
the random initial data y®(-) and y„(-) determined by their Fourier coefficients = Pi jdj and 
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bj = 0% jbj (see (fl~3")) for the explicit relation between the Fourier coefficients and the initial data), 
that is, 

+00 

y„(t, x) = Yl {^ a ^ X]t + PIM'^) M*)- (32) 
3=1 

This new constant G<r2nd(Xw) is called randomized observability constant. 

Similarly, making a random selection of all possible initial data for the Schrodinger equation 

T 

dV g dt I , (33) 



© we replace C^(Xu>) by 



(A^)e^(C) \ Jo 

££?Ajl«*l a =l 



+00 



X>; nv A 

3=1 



where (/3J) j6 ]n» denotes a sequence of independent Bernoulli random variables on a probability 
space (X,A,F). This corresponds to considering the randomized observability inequality 

c¥Hx«)\\y°\\ 2 D(A) < E ( / T / |9 t ^M| 2 dVWt] > (34) 



'0 v'CJ 



for every y°(-) £ -D(^4), where y„ denotes the solution of the Schrodinger equation with the ran- 
dom initial data y°(-) determined by its Fourier coefficients Cj — (3fcj (see (|22|) for the explicit 
dependence between the Fourier coefficients and the initial data), that is, 

+00 

y v {t,x)=Y J ^Oje iX ^<i> j {x). 
3=1 

The following theorem, whose proof is done in Section 12.41 provides one more motivation of 
studying the second problem (1111) . 

Theorem 2. There holds 

2 4^L(X") = C { T % d ( X .) = T inf f ^(x) 2 dV g = TJ( Xu ), 
for every measurable subset lu of Q. 

Remark 5. It can be easily checked that Theorem [2] still holds true when considering, in the 
above randomization procedure, more general real random variables that are independent, have 
mean equal to 0, variance 1, and have a super exponential decay. We refer to |12[ 1151 \W\ for more 
details on these randomization issues. Bernoulli and Gaussian random variables satisfy such ap- 
propriate assumptions. As proved in [17], for all initial data (j/°, y 1 ) £ D(A 1 / 2 ) x X , the Bernoulli 
randomization keeps constant the D^A 1 / 2 ) x X norm, whereas the Gaussian randomization gen- 
erates a dense subset of D(A 1 ' 2 ) x X through the mapping i?( y o i2/ i) : v G X h-> (y°,yj) provided 
that all Fourier coefficients of (y°, y 1 ) are nonzero and that the measure 9 charges all open sets of 
IR. The measure /z^o^i) defined as the image of V by R( y a y i^ strongly depends both on the choice 
of the random variables and on the choice of the initial data (y°, y 1 ). Properties of these measures 
are established in [T7] . 

Remark 6. It is easy to sec that C^Ld(Xw) ^ C^ r> (xJ) and C^ tallA (xu) > Ct\xu), for every 
measurable subset u of f2, and every T > 0. 
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Remark 7. As mentioned previously, the problem of maximizing the deterministic (classical) ob- 
servability constants Cfpfcw) and Crp(xw) defined by © and (J7J) respectively, over all possible 
measurable subsets u of fl of measure V g (oj) = LV g (Cl), is open and is probably very difficult. 
It can however be noticed that, for practical issues, it is actually more natural to consider the 
problem of maximizing the randomized observability constants defined by (I30[) and f|33|) respec- 
tively. Indeed, when considering for instance the practical problem of locating sensors in an optimal 
way, the optimality should be thought in terms of an average with respect to a large number of 
experiments. From this point of view, the deterministic observability constants are expected to 
be pessimistic with respect to their randomized versions. Indeed, in general it is expected that 

CZIM > C^\ Xu ) and 41nd0fe) > C { T S) {xu). 

In dimension one, with 51 = [0,7r] and Dirichlet boundary conditions, it follows from [54, Proposi- 
tion 2] (where this one-dimensional case is studied in detail) that these strict inequalities hold if 
and only if T is not an integer multiple of 7r (note that if T is a multiple of 2n then the equalities 
follow immediately from Parseval's Theorem). Note that, in the one-dimensional case, the GCC 
is satisfied for every T ^ 2tt, and the fact that the deterministic and the randomized observability 
constants do not coincide is due to crossed Fourier modes in the deterministic case. 

In dimension greater than one, there is a further class of obvious examples where the strict 
inequality holds. This is particularly the case when one is able to assert that C^ / \xu) = 
whereas J(xuj) > 0. Such examples have been given and discussed in Remark Note however 
that, in this multi-dimensional case, this is due to the fact that the GCC can fail but the minimal 
spectral trace over w is positive. This fact was also observed in |45j when characterizing the decay 
rates for multi-dimensional dissipative wave equations. 



2.3 Proofs of Theorem Q] and Corollary Q] 



We prove Theorem Q] and Corollary [T] only for \xu) (wave equation). The proof for C^^x^) 
(Schrodinger equation) follows the same lines. For the convenience of the reader, we first prove 
Theorem[T]in the particular case where all the eigenvalues of A g are simple (it corresponds exactly 
to the proof of Corollary [J) and we then comment the generalization to the case of multiple 
eigenvalues. 

From (fl"2"l) . we have y{t,x) — Yl~j=i yj{pi x ) with 

Vj (t, x) = ( aj e iX ^ + bje-^^ix). (35) 
Without loss of generality, we consider initial data (y ,?/ 1 ) € D(A 1 / 2 ) xX such that || (y°, y 1 )\\ 2 D ( A i/2' )x x 

2, in other words such that J2jew* -^(Kl 2 + \ b j\ 2 ) = 1 (using ([iUO . 
Setting 

^T[a,b) = — = — G T (x u ), 



T \\(y°,y 1 )\\ D{A i/2 )> 



X 



2T 



we write for an arbitrary N € IN* , 

r-T 



N 



+00 



N 



k=N+l 



+2Ke[^ % (i,x) yk(t,x)\ \ dV g dt. 

,j=l fe=Ar+l 



(36) 



Using the assumption that the spectrum of A consists of simple eigenvalues, we have the following 
result. 
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Lemma 1. With the notations above, 



lim — 



Juj 



N 
3 = 1 



N 

dv g dt = Y, \ 2 (KI 2 + N 2 ) / fa¥) 2 dv g 

3=1 Ju 



Proof. Since the sum is finite we can invert the infimum (which is a minimum) and the limit. Now, 
we write 

2 

TJo 



N 
3=1 



dV 9 dt = ^Y, X >33 l M*) 2 dV g 
3 = 1 Ju} 
1 N N 

+ f^2^2 X 3 X kajk / <t>j{x)(j>k{x) dVc 

A — 1 J. 1 « W 



J=l fc=l 



where is defined by p^|) . Using (fTT)) and (fT5j) . we get 



lim ^ = | a ,l 2 + |&,| 2 , 



for every j £ IN* and, using that the spectrum of A consists of simple eigenvalues, 

4maxi^^jv(A 3 , A fe ) 



\Otjk\ < 

whenever j ^ k. The conclusion follows easily. 
Let us now estimate the remaining terms 

l-T 



R 



1 
T 



3=N+1 



dV g dt 



(rp -\-QO \ 

f [ £>M J2 Vk(t,x)dV g dt) 
J ° Ju j=l k=N+l J 



and 



of the right-hand side of (1361) . 

Estimate of R. Using the fact that the 4>j's form a hilbertian basis, we get 

2 



R < 



1 

TJ 



+ 00 



=JV+1 



dV g dt 



e~ lX ^\ 2 dt 



j=N+l ■ 



i E AMT(Kf + |^| 2 )-i-5fte(aA 



g2iAjT _ j 



and finally 



i? 



3=N+1 



(37) 
□ 



(38) 



1G 



Estimate of 5. Using (fTT)) and the fact that Xj ^ Xk for every j e {1, , } and every 
k ^ N + 1 , we have 

with 



N 



S. 



N 



N +oo 

E E A ; - A; ' 

j=l k=N+l ■> 



T 



aj a k e i( - x >- Xk ^ sin (Xj - X k )- / ^(^(a:) dV £ 



AT +00 



XjXk 
1^ A~TaT j 

j=l k=N+l 3 K 



N +00 . . 



j=l k=N+l 

N +00 



A, + At 



< j a k e- i ^ +x "^ sin f (Aj + A fc )|) j ^(x)(j> k {x) dV g 



A A / T\ f 

E E -^6Ae- i(Aj - Afc)? sin (A. -Afe)- / ^{x)^ h {x)dV a 



Let us estimate S± . We write 

N „ +co 



iV 



E V, / E T^-e^-^T sin ( (A, - A fc )^ ) dF g 



fc=jv+i A -> ~ Afe 



and, using the Cauchy-Schwarz inequality and the fact that the integral of a nonnegative function 
over lo is lower than the integral of the same function over f2, one gets 



N 



X k a k 



+00 

y 



1/2 



E A W E T^S^foj-A*)? 



(A, - A fe )2 



1/2 



The last equality is established by expanding the square of the sum inside the integral, and by 
using the fact that the 4>kS are orthonormal in L 2 (Vl). Since the spectrum of A consists of simple 
eigenvalues (assumed to form an increasing sequence), we infer that Xk — Xj ^ Ajv + i — Ajv for all 
j e {!,••• ,N} and k > N + 1, and since A^jaj| 2 < 1, it follows that 



1 N / +oc 

\n+i - X N £j Vfc^+i 



1/2 



Afc|a fe | 



iV 



Aw+i — Ajv 



The same arguments lead to the estimates 



N 
Ajv ' 



N 
Xn ' 



N 



Xn+i — Xn ' 
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and therefore, 

1 \An Aat + i — Ajv, 



Now, combining Lemma [JJ with the estimates (|38|) and (p9|) yields that for every e > 0, there 
exist N e G IN* and T(e, N e ) > such that, if N ^ N £ and T ^ T(e, iV e ), then 



MM)-I>i(N 2 -H 6 il a ) / fcO*) 2 ^ 
j=i Jw 



As an immediate consequence, and using the obvious fact that, for every rj > 0, there exists 



N v G N* such that, if TV ^ N n then 



one deduces that 



£A 2 (hf + |6,| 2 ) / cj )o {xfdV g -Y,^{\ a 1 \ 2 + \ b 1 \ 2 ) / ^W 2 ^ 



lim E T (a,6) = ^A J 2 (| aj | 2 + |6 J f) / ^(xf dV g 



At this step, we have proved the following lemma, which improves the statement of Lemma [TJ 



Lemma 2. Denoting by aj and bj the Fourier coefficients of (yP^y 1 ) defined by (| 1 3[) . there holds 
lim IT/ |y(i,x)| 2 dK,di = ^A 2 (K| 2 + |6 J f) / ^{xf dV g . 
Corollary [1] follows, noting that 

+oo „ „ 

X>i(|aj| 2 + I*il 2 ) / ^(^) 2 rfV 9 = inf / ^-(s) 2 ^. 



+00 

inf 

££?A?(M a +IM 2 )=i 



To finish the proof, we now explain how the arguments above can be generalized to the case of 
multiple eigenvalues. In particular, the statement of Lemma 1 is adapted in the following way. 

Lemma 3. Using the previous notations, one has 

2 



lim — 



xeu 



/ 




2 




2 \ 




2J A k a k (j) k (x) 


+ 






\ 


fce/(A) 




kei(x) 


/ 



dV a 



Proof. Following the proof of Lemma [TJ simple computations show that 

2 



1 



N 



dVgdt = E AjAkUjk ( 4>j{x)4> k {x)dV g 

xeu 0',fe)e/(A) a 

+ ^ E E X i X kOtjk I <f>j{x)(j) k (x)dV g , 



A^/i k£l(p) 
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where 



li m a i k - [ a J ak + h i hk if C?' fc ) e J ( A ) 2 ' 
t->+oo T \ if j £ /(A), k £ with (A,/i) £ t/ 2 and A 7^ /i. 

The conclusion of the lemma follows. 



□ 



To derive Theorem [TJ it suffices to note that the previous estimates on R and 5 are still valid 
and that 



inf y 



y A fe a fe fc (a;) 



fceJ(A) 



inf 

(cf=)j e] N*e« 2 (C) 



E 

AGC/ 



c k (j) k (x) 

kei(x) 



dVa 



2.4 Proof of Theorem d 

From Fubini's theorem, using the fact that the random laws are independent, of zero mean and of 
variance 1, we have 



C { t%M = 2 inf 

(\ 2 Cj )e£ 2 (C) Jo 

£+«A}|«y| a =l 



+00 



dV g dt 



+OO 



^AfcyV^) 

3=1 

,fc=i 
. +00 

T inf J2\ c ^^) 2 dV g 



E£5|c,| a =i 



E^IcjI^i 
T inf / 4>AxfdV q 



The proof for C<j^2n d (Xw) i s similar. 



3 First problem: optimal design for fixed initial data 

This section is devoted to solving the first problem, that is, the problem of best observation for fixed 
initial data. Throughout the section, we fix initial data (y ,?/ 1 ) £ D(A 1 / 2 ) x X (resp., y° £ D{A)) 
for the wave equation (TTJ) (resp., for the Schrodinger equation fl2J)), and we consider their associated 
coefficients a^ , (i, j) £ (IN* ) 2 , defined by (TTB)) (resp., by ([25)1 ). The next considerations are valuable 
for both wave and Schrodinger equations. For every x £ f2, we define 

<p(x)= / |9 t t/(t,a;)| 2 rft = XiXjaij^i(x)4>j(x), (40) 
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where y is defined by (|12[) or ([^T|) and is such that y G C°(0, T; Z^A 1 / 2 )), and where the coefficients 
a>ij are defined by p^|) or (l2"51) . Note that the function (/? is integrable on ft. Then, from (|T5j) or 
(glD, there holds 

Gt(x^)= / <p(x)W g , (41) 

for every measurable subset u> of f2. 
3.1 Main result 

Theorem 3. There exists at least one measurable subset uj of £1, solution of the first problem, 
characterized as follows. There exists a real number A such that every optimal set uj is contained 
in the level set {tp ^ A}, where the function ip defined by (|40p is integrable on f2. 

Moreover, if M is an analytic Riemannian manifold, if f2 has a nontrivial boundary of class 
C°° and if there exists R > such that 

£ % (ll^Vlli- + \\A^V\\h) 1/2 < +oo, (42) 
in the case of the wave equation, and 

(43) 

3=0 J ' 

in the case of the Schrddinger equation, then the first problem has a uniqu^ solution Xu> where u> 
is a measurable subset of VL of measure LV g (Q), satisfying moreover the following properties: 

• to is semi-analyti^, and has a finite number of connected components; 

• if M — Ft" , if fl is symmetric with respect to an hyperplane, if y° o a — y° and y 1 o a = y 1 
where a denotes the symmetry operator with respect to this hyperplane, then uj enjoys the 
same symmetry property; 

• for Dirichlet boundary conditions, there exists n > such that d(uj,dtt) > n, where d denotes 
the Riemannian distance on M . 

The first statement of this theorem covers the case where <9f2 = 0. In this case, fi is a compact 
connected analytic Riemannian manifold. 

Proof of Theorem [3J The existence and the characterization in function of the level sets of ip of a set 
ui maximizing (|41l) is obvious since the function ip is integrable on Q. Let us prove the second part 
of the theorem, for the wave equation. First of all we claim that, under the additional assumption 

3 Similarly to the definition of elements of L p -spaces, the subset ui is unique within the class of all measurable 
subsets of f! quotiented by the set of all measurable subsets of f2 of zero measure. 

4 A subset u of a real analytic finite dimensional manifold M is said to be semi-analytic if it can be written in 
terms of equalities and inequalities of analytic functions, that is, for every x £ u), there exists a neighborhood U of 
x in M and 2pq analytic functions gij, hij (with 1 ^ i ^ p and 1 ^ j ^ q) such that 

v 

LOf\U = {J{y & U \ gij (y) = and h i3 {y) > 0, j = l,...,q}. 

i = l 

We recall that such semi-analytic (and more generally, subanalytic) subsets enjoy nice properties, for instance they 
are stratifiable in the sense of Whitney (see 1271 1341 ) . 



20 



(HH), the corresponding solution y of the wave equation is analytic over H + x 57. Indeed, we hrst 
note that the quantity 

\\Ai^y(t,-)\\l 2 + \\A^'*d t y(t,-)\\l 2 

is constant with respect to t. Then, since 57 has a smooth boundary, it follows from (l42l and from 
the Sobolev imbedding theorems that there exists C > such that 



\\ v W( t .)\\ ^C (2n + fc)! 
\\y [i, )\\oo 5* ^ ^ 2 «+fe > 

for every t ^ and every integer fc. The claim follows. As a consequence, the function ip defined 
by (|40p is analytic on 57. Hence </? cannot be constant on a subset of positive measure (otherwise 
by analyticity it would be constant on 57 and hence equal to due to the boundary conditions). 
This ensures the uniqueness of the optimal set u. 

The first additional property follows from the analyticity properties. The symmetry property 
(if M = R n ) follows from the fact that tp o a(x) = tp(x) for every x £ 57. If u) were not symmetric 
with respect to this hyperplane, the uniqueness of the solution of the first problem would fail, 
which is a contradiction. For Dirichlet boundary conditions, since tp(x) = f Q T \dty(t, x)\ 2 dt — for 
every x € 357, it follows that tp reaches its global minimum on the boundary of 57. □ 

Remark 8. The solution of the first problem depends on the initial data under consideration. 
More specifically, it depends on their associated coefficients , defined by (fT6|) in the case of the 
wave equation, and by (|25|) in the case of the Schrodinger equation. Note that there exist an infinite 
number of initial data (z°, z 1 ) € D{A 1 / 2 ) x X (resp. z° G D(A)) having the same coefficients ctij 
than (y ,?/ 1 ) (resp., y°), and all of them lead to the same solution of the first problem. Similar 
considerations have been discussed in the one-dimensional case in [53] . 

Remark 9. We have seen that the optimal solution may not be unique whenever the function tp 
is constant on some subset of 57 of positive measure. More precisely, assume that tp is constant, 
equal to c, on some subset I of 57 of positive measure |/|. If \{<p > c}\ < LV g (fl) < \{ip > c}\ 
then there exists an infinite number of measurable subsets w of 57 maximizing (|4ip , all of them 
containing the subset {tp > c}. The part of co lying in {tp = c} can indeed be chosen arbitrarily. 

Note that there is no simple characterization of all initial data for which this non-uniqueness 
phenomenon occurs, however to get convinced that this may indeed happen it is convenient to 
consider the one-dimensional case where T is moreover an integer multiple of 2tt (see Remark [3]), 
with Dirichlet boundary conditions. Indeed in that case the functional Gt does not involve any 
crossed terms and therefore the corresponding function tp reduces to tp{x) — ^j a jj snl2 (.? a; ); 

with atjj given by (TPS)) . Writing sin 2 (jo;) = | — ^cos(2jx) permits to write tp as a Fourier series 
whose sine Fourier coefficients vanish and cosine coefficients are nonpositive and summable (because 
of (|19p ). Hence, to provide an explicit example where the non-uniqueness phenomenon occurs, it 
suffices to consider a nonpositive triangle function defined on [? — a, 5 +a], for some a > 0, equal 
to outside. Its Fourier coefficients are the values on integers of the Fourier transform of the 
nonpositive triangle function, hence are negative and summable. The rest of the construction is 
obvious. We refer to [S3] for details on this one-dimensional case. 

Note that it is easy to generalize such a characterization in a n-dimensional hypercube for the 
Schrodinger equation, since in this case the solution remains periodic with period 27r and Gt does 
not involve any crossed terms. 

3.2 Further comments on the complexity of the optimal set 

It is interesting to raise the question of the complexity of the optimal sets solutions of the first 
problem. In Theorem [3] we prove that, if the initial data belong to some analyticity spaces, then 
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the (unique) optimal set u> is the union of a finite number of connected components. Hence, 
analyticity implies finiteness and it is interesting to wonder whether this property still holds true 
for less regular initial data. 

In what follows we show that, in some particular cases (those mentioned above where the 
periodicity of solutions of the wave and Schrodinger equations can be exploited), there exist C°° 
initial data for which the optimal set u has a fractal structure and, more precisely, is of Cantor 
type. 

Proposition 1. For the one- dimensional wave equation on [0, 7r] (resp. for the n- dimensional 
Schrodinger equation on the hypercube [0,Tr] n ) with Dirichlet boundary conditions, there exist C°° 
initial data (y° ,y l ) (resp. C°° initial data y° ) for which the first problem has a unique solution uj 
with fractal structure and thus, in particular, it has an infinite number of connected components. 

The proof of this proposition is quite technical and relics on a careful Fourier analysis construc- 
tion. It is done in Appendix lAl 



3.3 Several numerical simulations 

We provide hereafter a numerical illustration of the results presented in this section. 

According to Theorem [3l the optimal domain is characterized as a level set of the function ip. 
Some numerical simulations are provided on Figure [TJ with Q = [0,7r] 2 , L — 0.6, T = 3, y 1 = 
and 

N 

y°(x)— a n , k sm(nxi) sm(kx 2 ), 

n,k—l 

where Nq G IN* and (a n ^)n,ke'N* are real numbers. The level set is numerically computed using a 
simple dichotomy procedure. 

4 Second problem: uniform optimal design 

4.1 Preliminary remarks 

We define the set 

Ul = {Xuj | w is a measurable subset of f2 of measure V g (w) = LV g (£l)}. (44) 
Recall that the second problem (JTTJ) is written as 

sup J{Xw), (45) 

with 



The criterion J(Xui) can be seen as a spectral energy concentration criterion. For every j € IN*, 
the integral J <fij(x) 2 dV g is the energy of the j th eigenfunction restricted to to, and the problem is 
to maximize the infimum over j of these energies, over all subsets u> of measure V g {uS) = LV g (il). 

Since the set Ul does not have compactness properties ensuring the existence of a solution of 
(|45|) . we consider the convex closure of Ul for the weak star topology of L°° , 

U L = ja e L°°(ft, [0, 1]) | j a(x) dV g = LV g {Sl)X . (46) 
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Some level sets of y° 



Optimal domain for L=0.6 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



Figure 1: On this figure, Q — [0,7r] 2 with Dirichlet boundary conditions, L = 0.6, T = 3 and 
y = 0. At the top: No = 15 and a m ,t = n -±^ k -i ■ At the bottom: Nq = 15 and a n ,h = — — ■ 
On the left: some level sets of y°. On the right: representation of the optimal domain for the 
corresponding choice of y . 



This convexification procedure is standard in shape optimization problems where an optimum may 
fail to exist because of hard constraints (see e.g. [TU]). 

Replacing G Wi with a g Ul, we define a convexified formulation of the second problem 
(g5|) by 

sup J (a), (47) 
aeu L 

where 

J(a)= inf / a{x)4>Ax) 2 dV g . (48) 
Jo 

Since J(a) is defined as the infimum of linear continuous functionals for the weak star topology of 
L°°, it is upper semi continuous for this topology. This yields to the following result. 



Lemma 4. The problem (I47|) has at least one solution. 
Obviously, there holds 

sup inf / ^(i)^, ,(x) 2 dV g < sup inf / a(x)<j)j(x) z dV„. (49) 
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Note that, since the constant function a(-) = L belongs to Ul, it follows that sup J(a) L. In 

aEU L 

the next section, under an additional ergodicity assumption, we compute the optimal value (|47l) 
of this convexified problem and investigate the question of knowing whether the above inequality 
is strict or not. In other words we investigate whether there is a gap or not between the problem 
(|4"5j) and its convexified version (|4"7|) . 

Remark 10. Comments on the choice of the topology. 

In our study we consider measurable subsets u> of il, and we endow the set L°°(il, {0, 1}) of all 
characteristic functions of measurable subsets with the weak-star topology. Other topologies are 
used in shape optimization problems, such as the Hausdorff topology. Note however that, although 
the Hausdorff topology shares nice compactness properties, it cannot be used in our study because 
of the measure constraint on u. Indeed, the Hausdorff convergence does not preserve measure, 
and the class of admissible domains is not closed for this topology. Topologies associated with 
convergence in the sense of characteristic functions or in the sense of compact sets (see for instance 
[321 Chapter 2]) do not guarantee easily the compactness of minimizing sequences of domains, 
unless one restricts the class of admissible domains, imposing for example some kind of uniform 
regularity. 

Remark 11. We stress that the question of the possible existence of a gap between the original 
problem and its convexified version is not obvious and cannot be handled with usual T-convergence 
tools, in particular because the function J defined by (pf5)l is is not lower semi-continuous for the 
weak star topology of L°° (it is however upper semi-continuous for that topology, as an infimum 
of linear functions). To illustrate this fact, consider the one-dimensional case of Remark [3l In this 

specific situation, since <j>j (x) = \l — sin(jx) for every j £ IN* , one has 



and the supremum is reached with the constant function a(-) = L. Consider the sequence of subsets 
wjv of [0,7r] of measure Lit defined by 



for every N € IN*. Clearly, the sequence of functions \ui N converges to the constant function 
a(-) = L for the weak star topology of L°°, but nevertheless, an easy computation shows that 




for every a € Ul- Since the functions x i-> sin 2 (ja;) converge weakly to 1/2, it clearly follows that 
J (a) ^ L for every a £ Ul- Therefore, 



sup J(a) = L, 

aGU L 




f 



sin 2 (jx) dx 




-|sinKf 1 if (N + l) \j, 
+ j7 sin ( ijf- ) otherwise, 



and hence, 
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This simple example illustrates the difficulty in understanding the limiting behavior of the func- 
tional because of the lack of the lower semicontinuity, what makes possible the occurrence of a gap 
in the convexification procedure. In Section 14.21 we will prove that there is no such a gap under 
an additional geometric spectral assumption. 



4.2 Main results 

In what follows, we make the following assumptions on the basis (0|)j 6 ]N» of eigenfunctions under 
consideration. 

Weak Quantum Ergodicity on the basis (WQE) property. There exists a 
subsequence of the sequence of probability measures fij — (j) 2 dV g converging vaguely to 
the uniform measure v ^ dV g . 

Uniform L°° -boundedness property. There exists A > such that 

H0jllz<»(n) < A > (5°) 

for every j £ JV* . 

Note that the two assumptions above imply what we call the L°° - Weak Quantum Ergodicity on 
the base (L°°-WQE) property^, that is, there exists a subsequence of {<p~)j^' converging to y-jrr: 
for the weak star topology of L°°(il). 
Obviously, this property implies that 



sup inf / a(x)4>j(x) dV g = L, (51) 
aeu L J' e]N * Jn 

and moreover the supremum is reached with the constant function a — L on il. 

Remark 12. In general the convexified problem (|47|) does not admit a unique solution. Indeed, 
under symmetry assumptions on Q there exists an infinite number of solutions. For example, in 
dimension one, with il = [0, n], all solutions of (H71) are given by all functions oUAl whose Fourier 
expansion series is of the form a{x) = L + J2^=i( a j cos(2jx) + bj sin(2jx)) with coefficients aj ^ 0. 

It follows from (|49j) and ([51]) that 

sup inf / 6Ax) 2 dV a ^ L. 

The next result states that this inequality is actually an equality. 

Theorem 4. // the WQE and uniform L°° -boundedness properties hold, then 



sup inf / 4>j(x) 2 dV = L, 



(52) 



for every L £ (0, 1). In other words, under these assumptions there is no gap between the original 
problem (|45l) and the convexified one. 



5 The wording used here is motivated and explained further in a series of remarks. 
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It follows from this result, from Corollary [T] and Theorem [51 that the maximal value of the 
randomized observability constants 2C^2 nd (xw) = C-Trandix^) over the set Ul is equal to TL, 
and that, if the spectrum of A is simple, the maximal value of the time asymptotic observability 

constants 2C^\x u ) — Coo\xw) over tne set Ml 1S equal to L. 

The question of knowing whether the supremum in (|52p is reached (existence of an optimal set) 
is investigated in Section I4T41 

Theorem [4] is established within the class of measurable subsets. We next state a similar (but 
distinct) result within the class of measurable subsets whose boundary is of measure zero. We 
define the set 

U h L = {xui G U L | Vg(duj) = 0}. (53) 

This is the set of all characteristic functions of Jordan measurable subsets of of measure LV g {Q). 
We make the following assumptions. 

Quantum Unique Ergodicity on the base (QUE) property. The whole sequence 
of probability measures Uj = (f>j dV g converges vaguely to the uniform measure v ^ dV g . 

Uniform L p -boundedness property. There exist p £ (1, +oo] and A > such that 

Hjh'pw <A (54) 

for every j € JV* . 

Theorem 5. Assume that dfl is Lipschitz whenever it is nonempty. If the QUE and uniform 
LP -boundedness properties hold, then 

sup inf / (j)j(x) 2 dV g = L, (55) 

for every L € (0, 1). 

Theorems |4] and [5] are proved in Sections 14.51 and 14.61 respectively. 

Remark 13. It follows from the proof of Theorem[5]that this statement holds true as well whenever 
the set U\ is replaced with the set of all measurable subsets w of fi, of measure V g (ui) = LV g (£l), 
that are moreover either open with a Lipschitz boundary, or open with a bounded perimeter. 

Remark 14. The assumptions made in Theorems 2] or [5] are sufficient conditions implying (I5"2"j) 
or (|55p. but they are however not sharp, as proved in the next proposition. 

Proposition 2. Consider the Dirichlet-Laplacian on the domain Q defined as the unit disk of the 
Euclidean two-dimensional space. Then, for every p £ (l,+oo] and for any basis of eigenf unctions, 
the uniform L p -boundedness property is not satisfied, and QUE does not hold as well. However, 
the equalities (|52[) and (|55j) hold true. 

To establish this result, in the proof of this proposition (done in Section |4~7]) we use the explicit 
expression of certain semi-classical measures in the disk (weak limits of the probability measures 
(pj dV g ). Among these quantum limits, one can find the Dirac measure along the boundary which 
causes the well known phenomenon of whispering galleries. Having in mind this phenomenon, 
it might be expected that there exists an optimal set, concentrating around the boundary. The 
calculations show that it is not the case, and (f52|) and (f55|) are proved to hold. 

The next section is devoted to gather some comments on the ergodicity assumptions made in 
these theorems. 
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4.3 Comments on ergodicity assumptions 

This section is organized as a series of remarks. 

Remark 15. The assumptions of Theorem 0] hold true in dimension one. Indeed, it has already 
been mentioned that the eigenfunctions of the Dirichlet-Laplacian operator on Cl = [Q,w] are given 

by 4>j(x) = y^sba.{jx), for every j 6 IN*. Therefore clearly the whole sequence (not only a 

subsequence) (0|)j- e ]N* converges weakly to ^ for the weak star topology of L°°(0,7r). The same 
property clearly holds for all other boundary conditions considered in this article. 

Remark 16. In dimension greater than one the situation is more intricate, but we have the 
following facts. 

Any hypercube (tensorised version of the previous one-dimensional case) or flat torus satisfies 
the assumptions. Indeed, the whole sequence of eigenfunctions is uniformly bounded and converges 
to a constant for the weak star topology of L°° . 

Generally speaking, these assumptions are related to ergodicity properties of Cl. Before provid- 
ing precise results, we recall the following well known definition. 



Quantum Ergodicity on the base (QE) property. There exists a subsequence of 
the sequence of probability measures fXj = (fq dV g of density one converging vaguely to 
the uniform measure v ^ dV g . 

Here, density one means that there exists T C IN* such that 

hm *iiHli^Il = 1 . 

Obviously, QE implies WQE@. The well known Shnirelman Theorem asserts that the property QE is 
satisfied on every compact ergodic Riemannian manifolds having no boundary (see [2"T1 1551 [5^1 IBT] ) . 
For domains having a boundary, it is proved in [55] that, if the domain Cl is a convex ergodic billiard 
with W 2 '°° boundary, then the property QE is satisfied. The result of [26] has been extended to 
arbitrary ergodic manifolds with piecewise smooth boundaries in |71] and [29] . 

Note that this result lets however open the possibility of having an exceptional subsequence 
of measures fj,j converging vaguely to something else. We will come back on this interesting issue 
later. 

Actually these results relating the ergodicity of CI (seen as a billiard where the geodesic flow 
moves at unit speed and bounces at the boundary according to the Geometric Optics laws) to 
the QE property are even stronger, for two reasons. Firstly, they are valid for any Hilbertian 
basis of eigenfunctions of A, whereas here we make this kind of assumption only for the specific 
basis ((^j)jg]N* that has been fixed at the beginning of the article. Secondly, they establish that a 
stronger microlocal version of the QE property holds for pseudodifferential operators, in the unit 
cotangent bundle S*Cl of CI: more precisely it is proved that, under ergodicity assumptions on the 
manifold, a density one subsequence of the linear functionals Pj(P) = (P<f>j,(f)j), defined on the 
space of zero-th order pseudo-differential operators P, converges vaguely to the uniform Liouville 
measure. It is a much stronger conclusion since it says that the eigenfunctions become uniformly 
distributed on the phase space S*Cl and not just on the configuration space CI. Here however we 
do not need (de)concentration results in the full phase space, but only in the configuration space. 
This is why, following [BU] , we use the wording "on the base" . 



"Note that, up to our knowledge, the notion of WQE is new, whereas the notions of QE and QUE are classical 
in mathematical physics. 
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Note that the vague convergence of the measures (Xj is weaker than the convergence of the 
functions (j> 2 for the weak star topology of L°°(SY). Since il is bounded, the property of vague 
convergence is equivalent to saying that, for a subsequence of density one, J M <fij(x) 2 dV g converges 
to V g (ui)/V g (fl) for every measurable subset u> of £1 such that V g (duj) = (Portmanteau theorem). 
In contrast, the property of convergence for the weak star topology of L°°(fl) is equivalent to 
saying that, for a subsequence of density one, J <j>j{x) 2 dV g converges to V g (ui)/V g (O) for every 
measurable subset oj of Q. Under the assumption that all eigenfunctions are uniformly bounded 
in L°°(fi), both notions are equivalent. This is the case for instance in flat tori. But, for instance, 
if Q is a ball or a sphere of any dimension, then the eigenfunctions of the Laplacian are not 
uniformly bounded. This is well known to be a delicate issue (see [69]). We refer to [70] where 
it is conjectured that flat tori are the sole compact manifolds without boundary where the whole 
family of eigenfunctions is uniformly bounded in L°° . 

Note that the notion of L°°-QE property, meaning that the above QE property holds for the 
weak star topology of L°°, is defined and mentioned in [5j5] as a delicate open problem. As 
said above we stress that, under the assumption that all eigenfunctions are uniformly bounded in 
L°°(fi), QE and L°°-QE are equivalent. 

To the best of our knowledge, nothing seems to be known on the uniform L p -boundedness 
property (|50p . This property holds for fiat tori but does not hold for balls or spheres. 

Remark 17. Let us comment on the QUE property, which is an important issue in quantum and 
mathematical physics. Note indeed that the quantity J u (j> 2 (x) dV g is interpreted as the probability 
of finding the quantum state of energy A 2 in u>. We stress again on the fact that, here, we consider a 
version of QUE in the configuration space only, not in the full phase space. Moreover, we consider 
the QUE property for the basis (<f>j)jevi* under consideration, but not necessarily for any such 
basis of eigenfunctions. 

First of all, QUE obviously holds true in the one-dimensional case of Remark[3] (see also Remark 
Hip but it does however not hold true for multi-dimensional hypercubes. 

In the general multi-dimensional case, many interesting open questions and issues occur. As 
in Remark I16( consider for every j G IN* the probability measure fj,j = (fq dV g , representing in 
quantum mechanics the probability of being in the state tj>j (or, probability density of finding a 
particle of energy A 2 at x). An interesting question is to know whether the supports of these 
measures tend to equidistribute or can concentrate as j — > +oo. As already mentioned, under 
ergodicity assumptions, the property QE holds true, that is, a subsequence of density one of 
(MjOjeW converges vaguely to the uniform measure on Q. But this result lets open the possibility 
of having an exceptional subsequence converging to some other measure. Typically it may happen 
that a subsequence of density zero converges to an invariant measure like for instance a measure 
carried by closed geodesies. These so-called (strong) scars result from such energy concentration 
phenomena, that are allowed in the context of Shnirelman Theorem. This fascinating question of 
knowing whether the quantum states of such an ergodic system can concentrate or not on some 
instable closed orbits or on some invariant tori generated by such geodesies is still widely open in 
mathematics and physics. We refer to [9, 25J for results showing a scarring phenomenon on the 
periodic orbits of the dynamics of the quantum Arnold's cat map. Note however that, as already 
mentioned, here we arc concerned with concentration results in the configuration space only. 

The QUE property on the base, stating that the whole sequence of measures fij = <fi 2 dV g con- 
verges vaguely to the uniform measure, postulates that there is no such concentration phenomenon 
(see [55]). 

We recall that what is called the billiard in f2 is the dynamical system posed on the unit 
cotangent bundle of f2, representing (almost all) trajectories in il along geodesies with unit speed 
and reflecting on the boundary of fi according to the usual reflection rules. The ergodicity of fi is 
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Figure 2: Top left: a Barnctt billiard, conjectured in [6] to satisfy QUE. Top right: a Bunimovich 
stadium satisfying QE but not QUE for almost every value of the straight edge t. Bottom left 
and right: two shapes with piecewise regular boundary that are not ergodic and thus not quantum 
ergodic. 



however just a necessary assumption for QUE to hold. Note that strictly convex billiards whose 
boundary is C 6 are not ergodic in the phase space (see [43]), and there are sequences of positive 
density of eigenfunctions which concentrate on caustics. It has been shown in |38) that rational 
polygonal billiards are not ergodic in the phase space, while polygonal billiards are generically 
ergodic. In [33] the authors prove that quantum ergodicity on the base holds in any rational 
polygon^. It has been proved recently in [28] that there exist some convex sets satisfying QE 
but not QUE (independently on the basis of eigenfunctions under consideration). More precisely, 
in this reference the author studies the particular case of a stadium St with straight edge t (see 
Figure [2] top right). He shows that for almost every t the stadium is not quantum unique ergodic 
although it is quantum ergodic. He exhibits some particular quantum limits giving a positive 
mass on the set of bouncing ball trajectories. Up to now obtaining sufficient conditions on the 
domain such that QUE holds is a widely open difficult question. It was conjectured in [56 that 
every compact negatively curved manifold satisfies QUE (see |57] for a recent survey) . A numerical 
method has been developed in [6] in order to compute the 700000 first modes for a (planar) domain 
analogue of variable negative curvature, and the results confirm numerically the QUE conjecture 
for such general systems (see Figure [2] top left). Note that the QUE property has been proved to 
hold on arithmetic manifolds in 46 . Finally, note that, using a concept of entropy, the authors of 
[2j |3] show that on a compact manifold of negative curvature the eigenfunctions cannot concentrate 
entirely on closed geodesies and at least half of their energy remains chaotic (see also [18] for related 

7 A rational polygon is a planar polygon whose interior is connected and simply connected and whose vertex 
angles are rational multiples of w. 
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issues). Up to now, except the case of arithmetic manifolds, there does not exist any example of 
multi-dimensional domain in which QUE holds, and this is still currently one of the deepest issues 
in mathematical physics. 



Finally, having Figure [3] in mind , it is not surprising that the stable trapped ray of the left-side 
figure causes a concentration of eigenf unctions, due to the convexity of the domain f2. On the right- 
side figure, we have shaped a neighorhood of a domain Q in which negative curvature is suggested 
by the hyperbolic boundary, and there is a unique trapped ray, which is instable. Due to this 
instability feature, it might be expected that the energies of eigenfunctions spread away and that 
QUE holds true. This intuition is however not true. In [52], the authors build a compact surface 
of R 3 endowed with a metric of negative curvature, by truncating (for instance) an hyperboloid 
symmetrically with respect to its center, and considering Dirichlet boundary conditions on both 
truncated sides. Then, they show the existence of sequences of eigenfunctions of the Dirichlct- 
Laplacian that concentrate around the equator that is a closed instable geodesic. This surprising 
example shows that the QUE property, as well as the QUE conjecture, is definitely a global one, 
and cannot be inferred from local considerations. 

It can however be noticed that, as a consequence of [2J, the arc-length measure along a closed 
geodesic on a negatively curved manifold cannot be a quantum limit (see also |57)V 

Remark 18. The results Theorems [4] and [5] are similar but distinct. The QUE property assumed 
in Theorem [5] is a very strong one and as said above up to now examples of domains in dimension 
more than one satisfying QUE are not known. The proofs of these results, provided in Sections 14. 51 
and 14. 61 are of a completely different nature. In particular, our proof of Theorem U is short but 
does not permit to get an insight on the possible theoretical construction of a maximizing sequence 
of subsets. In contrast, our proof of Theorem [5] is constructive and provides a theoretical way of 
building a maximizing sequence of subsets, by implementing a kind of homogenization procedure. 
Moreover, this proof highlights the following interesting feature: 

It is possible to increase the values of J by considering subsets having an increasing 
number of connected components. 

Note that another way of building maximizing sequences is provided in Section [SJ by considering 
an appropriate spectral approximation of the problem, suitable for numerical simulations. 

Remark 19. The question of knowing whether there exists an example where there is a gap 
between the convexified problem (|47|) and the original one (J45J , is an open problem. We think 
that, if such an example exists, then the underlying geodesic flow ought to be completely integrable 
and have strong concentration properties. 




stable trapped ray 



instable trapped ray 



Figure 3: Stable and instable trapped rays 
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Note that, according to [37], the set of quantum limits on the unit sphere S 2 of H 3 is equal to 
the whole convex set of invariant probability measures for the geodesic flow that are time-reversal 
invariant (that is, invariant under the anti-symplectic involution (%, £) i— ► (a;, — £) on T*M). In 
particular, the Dirac measure along any great circle 7 on S 2 (defined as an equator, up to a rotation) 
is the projection of a semi-classical measure. However, as already mentioned in our framework we 
have fixed a given basis (<pj)jeN* of eigenvectors, and we consider only the weak limits of the 
measures tfij(x) 2 dV g , whereas the result of [37] holds when one considers the limits over all possible 
bases. With a fixed given basis, we are not aware of any example having concentration properties 
strong enough to derive a gap statement. We refer to Section [4.81 and in particular to Proposition 
[3] for an example of a gap for an intrinsic variant of the second problem. 

Note that the usual basis on S 2 , consisting of spherical harmonics, does not satisfy the QE 
property. But, due to the high multiplicity of eigenvalues, there is an infinite dimensional mani- 
fold of orthonormal bases of eigenf unctions. Indeed, using the spherical harmonics, the space of 
orthonormal bases is identified with the infinite product rifc^o U(2k + 1) of unitary groups, and 
thus inherits of the corresponding probability Haar measure. It is then proved in [6 8) that, on the 
standard sphere, almost every orthonormal basis of eigenfunctions of the Laplacian satisfies the 
QE property. 

Remark 20. Our results here show that shape optimization problems are intimately related 
with the ergodicity properties of fi. Notice that, in the early article [20], the authors suggested 
such connections. They analyzed the exponential decay of solutions of damped wave equations. 
Their results reflected that the quantum effects of bouncing balls or whispering galleries play an 
important role in the success of failure of the exponential decay property. At the end of the article, 
the authors conjectured that such considerations could be useful in the placement and design of 
actuators or sensors. Our results of this section provide precise results showing these connections 
and new perspectives on those intuitions. In our view they are the main contribution of our article, 
in the sense that they have pointed out the close relations existing between shape optimization 
and ergodicity, and provide new open problems and directions to domain optimization analysis. 

4.4 On the existence of an optimal set 

In this section we comment on the problem of knowing whether the supremum in (|52p is reached 
or not, in the framework of Theorem 2] This problem remains essentially open except in several 
particular cases. 

For the one-dimensional case already mentioned in Remarks [3J[TT] and [T5] we have the following 
result. 

Lemma 5. Assume that Q = [0,tt]. Let L £ (0, 1). The supremum of J overUh (which is equal to 
L) is reached if and only if L = 1/2. In that case, it is reached for all measurable subsets uj C [0, n] 
of measure tt/2 such that uj and its symmetric image ui' = n — u) are disjoint and complementary 
in [0, 7r] . 

Proof. Although the proof of that result can be found in |3D] and in [51] , we recall it here shortly 
since similar arguments will be used in the proof of the forthcoming Lemma [B] 

A subset uj C [0, n] of Lebesgue measure Lit is solution of (|52"]) if and only if J sin 2 (jx) dx ^ 
Liv/2 for every j g IN*, that is, / cos(2ja;) dx ^ 0. Therefore the Fourier series expansion of \^ 
on [0, 7r] must be of the form 

+00 

L + J^(aj cos(2jx) + bj sm(2jx)), 
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with coefficients aj ^ 0. Let uj' = it — u> be the symmetric set of uj with respect to 7r/2. The 
Fourier series expansion of Xw' is 

+00 

L + cos(2jx) — sin(2j:c)). 

j=i 

Set g(:r) — L — 7}(xuj(x) + Xu'{ x ))i f° r almost every a; € [0,7r]. The Fourier series expansion of 
g is — a j cos(2j:c), with Oj ^ for every j g IN*. Assume that L =/= 1/2. Then the sets 

uj and w' are not disjoint and complementary, and hence g is discontinuous. It then follows that 
Yl'jLi a j — ~~ 00 ■ Besides, the sum Y^jLi a j * s a ^ so the limit of X^fcS a kA n {k) asn-) +00, where 
A„ is the Fourier transform of the positive function A„ whose graph is the triangle joining the 
points (— i,0), (0, 2n) and (-^,0) (note that A„ is an approximation of the Dirac measure, with 
integral equal to 1). This raises a contradiction with the fact that 

+00 

g(t)A n (t)dt = ^2a k A n (k), 
fe=i 

derived from Plancherel's Theorem. □ 

For the two-dimensional square O = [0, it} 2 studied in Proposition[5]we are not able to provide a 
complete answer to the question of the existence. We are however able to characterize the existence 
of optimal sets that are a Cartesian product. 

Lemma 6. Assume that il — [0,7r] 2 . Let L € (0, 1). The supremum of J over the class of all 
possible subsets uj = u)\ x uji of Lebesgue measure Lit 2 , where ui\ and ui2 are measurable subsets of 
[0,7r] ; is reached if and only if L € {1/4,1/2,3/4}. In that case, it is reached for all such sets uj 
satisfying 

■^(Xu,(x, y) + x w (tt - x, y) + Xu(x, tt - y) + Xu,^ - x,ir - y)) = L, 
for almost all (x,y) £ [0,7r 2 ]. 

Proof. A subset uj C [0, 7r] 2 of Lebesgue measure Lit 2 is solution of (|52l) if and only if 

— ^ / sin 2 (jx) sin 2 (fey) dx dy ^ L, 

for all (j,k) € (IN*) 2 , that is, 

cos(2j x) cos(2ky) dx dy ^ / cos(2jx)dxdy+ / cos(2ky) dx dy. (56) 



Set £ x = J* Xu{ x >y) dy for almost every x G [0, tt], and £ y — J* Xui(x,y) dx for almost every 
y G [0, 7t]. Letting either j or k tend to +00 and using Fubini's theorem in (|56p leads to 



/>7T />7T 

/ l x cos(2jcc) dx ^ and / 
Jo Jo 



£ y cos(2ky) dy ^ 0, 



for every j G IN* and every k € INT* . 

Now, if uj — uj\ x uj 2 , where uj\ and uj 2 are measurable subsets of [0, tt], then the functions 
x 1 y £ x and y^£ y must be discontinuous. Using similar arguments as in the proof of Lemma[5j it 
follows that the functions x 1— > £ x + £-k-x and y 1— > £ y + l v - v must be constant on [0, it], and hence, 



/•7T /• 7T 

/ l x cos(2jx) dx = and / 
Jo Jo 



£ y cos(2ky) dy = 0, 
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for every j g IN* and every k g IN* . Using (|56p , it follows that 



cos(2ja:) cos(2fcy) dy ^ 0, 
for all (j,k) G (IN*) 2 . The function .F defined by 

F ( x i y) = \(Xu{x, v) + Xw(tt - x, y) + XwO, 7T - y) + \^ ~ x, tt - y)), 

for almost all (x,y) £ [0,7r] 2 , can only take the values 0, 1/4, 1/2, 3/4 and 1, and its Fourier series 
is of the form 

L H — - ( f cos( 2 j u) cos(2kv) dudv) cos( 2 jx) cos(2ky), 

and all Fourier coefficients are nonnegative. Using once again similar arguments as in the proof 
of Lemma [5] (Fourier transform and Plancherel's Theorem), it follows that F must necessarily be 
continuous on [0,7r] 2 and thus constant. The conclusion follows. □ 

Remark 21. All results of this section can obviously be generalized to multi-dimensional domains 
f2 written as N cartesian products of one-dimensional sets. 

Remark 22. According to Lemma El if L = 1/2 then there exists an infinite number of optimal 
sets. Four of them are drawn on Figure |U It is interesting to note that the optimal sets drawn 
on the left-side of the figure do not satisfy the Geometric Control Condition mentioned in Section 
11.11 and that in this configuration the (classical, deterministic) observability constants \xu) 
and CjP(x u ) are equal to 0, whereas, according to the previous results, there holds 

This fact is in accordance with Remarks 2] and [7] 

Remark 23. Similar considerations hold for the two-dimensional unit disk. Actually it easily 
follows from Lemma [5] and from the proof of Proposition [2] that, for L = 1/2, the supremum of J 
over Ul is reached for every subset ui of the form 

Lo = {(r,8)e [0,1] x [0,2tt] | 9 e lo 6 }, 

where uig is any subset of [0, 27r] such that uig and its symmetric image ui' g = 2tt — uig are disjoint 
and complementary in [0, 2~k\. But we do not know whether or not there are other maximizing 
subsets. 

Remark 24. In view of the results above one could expect that when 17 is the unit iV-dimensional 
hypercube, there exists a finite number of values of L 6 (0,1) such that the supremum in (|52[) is 
reached. The same result can probably be expected for generic domains f2. But these issues are 
open. 

4.5 Proof of Theorem [4] 

Since we already have the inequality 

sup inf / Xu{x)<l>j{x) 2 dV g < L, 
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it suffices to prove that, for every e > 0, there exists Xu S Ul such that 



(j>j(x) 2 dV g - L 



for every j £ IN*. To prove this fact, we consider the function / defined by f(x) = (<pj(x) 
for every x £ SI. Using the fact that the eigenfunctions are uniformly bounded in L°°(f2), it 
is clear that f(x) £ £°°, for every x £ fi. Then, clearly, / G L 1 (ri,^°°) (using the Bochner 
integral), and j Q f dV g is the constant sequence of £°° equal to 1. For every e > 0, there exists a 
partition of Jl = Uj? =1 f2fe, with fife measurable, such that J Q \\(f — f n )\\e°°dV g < £/{L + 1), with 
/n = 12k=i a kXfi k - For every k £ {1, ...,n.}, let be a measurable subset of fi^ such that 
Ig(wfc) = LV^(Ofc). We set a; = U^ =1 Wfe. Note that, by construction, one has Xu> & Ml, and 

~ n r. n 

/ {Xm - L)U dV g = Y>fe / { X u - L) X n k dV g = V a fc (V 9 (w fc ) - LV g (n k )) = 0. 
Jn ,. Jq 



k=l 



k=l 



Therefore, there holds 



fdV g -L / fdV g 



( X u-L)fdV g 



(Xw - L)fndV g 



(Xu-L)(f - f n )dV g 



and the conclusion follows. 
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4.6 Proof of Theorem 

In what follows, for every measurable subset to of f2, we set 

Ij(u>)= I <f> 3 (x) 2 dV g , 

J uj 

for every j £ IN* . By definition, there holds 

J(lo) — inf Ij(uj). 

Note that it follows from QUE and from the Portmanteau theorem (see Remark ITB]) that, for every 
measurable subset to of O such that V g (uj) = LV g (£l) and V g (duj) — 0, one has Ij(u) — > L as 
j — » +00, and hence J(w) ^ L. 

Let oj be an open connected subset of of measure LV g (ft) having a Lipschitz boundary. In 
the sequel we assume that J(ujq) < L, otherwise there is nothing to prove. Using QUE, there exists 
an integer jo such that 

I j (u )>L-~(L-J(oj )), (57) 

for every j > j Q . Our proof below consists of implementing a kind of homogenization procedure 
by constructing a sequence of open subsets u>k (starting from loq) such that V g (uik) = LV g (uik) and 
lim J(cjfe) = L. Denote by ZJ^ the closure of w , and by ujq the complement of u>q in Q. Since 

k— S- + 00 

fl and wo have a Lipschitz boundary, it follows that wo and £7\u;o satisfy a 5-cone propertjH, for 
some 5 > (see [32l Theorem 2.4.7]). Consider partitions of ZJo and wg, 

if if 
u = \jFi and wg = (J F h (58) 

i=l i=l 

to be chosen later. As a consequence of the <5-cone property, there exists c$ > and a choice of 
partition (F^i^^k (resp. (-F<)i<i<.K") sucn that, for V g (Fi) small enough, 

Vie {l,--- ,K} (resp. Vie{l, - , - m - > c s (resp. ^ „ > ej , (59) 

V / diamF, V diam.Fi / 

where 7^ (resp., rji) is the inradiu^|of Fi (resp., Fi), and diami^ (resp., diami^) the Riemannian 
diameter of Fi (resp., of Fi). 

It is then clear that, for every i E {1, . . . , K} (resp., for every i £ {1, . . . , K}), there exists € Fi 
(resp., £ € ^) such that 5(^,^/2) Cf.C B{^,rn/c s ) (resp., 5(^,^/2) c F< C B(£i,r)i/c s )), 
where the notation i?(£, r/) stands for the open Riemannian ball centered at £ with radius 77. These 
features characterize a substantial family of sets (also called nicely shrinking sets), as is well known 
in measure theory. By continuity, the points £j and £i are Lebesgue points of the functions <fi 2 , for 
every j ^ jq. This implies that, for every j $J jo, there holds 

4> 3 {x) 2 dV g = VgiFi)^) 2 + o(V g (Fi)) as m -> 0, 

F 4 



°We recall that an open subset £7 of TR n verifies a <5-cone property if, for every x S d£l, there exists a normalized 
vector gj, such that C(y,t;x,8) C f2 for every t/ESlfl B(x,6), where C(y,£ x ,S) = {z G R" | (z — j/, §) ^ cos5||z — 
and < \\z — y\\ < 8}. For manifolds, the definition is done accordingly in some charts, for 5 > small enough. 
9 In other words, the largest radius of Riemannian balls contained in Fi. 
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for every i E {1, . . . , K}, and 



_ fa(xY dV g = VgiF^j&y + o(V g (Fi)) as ft -> 0, 
for every i G {1, . . . , if}. Setting 77 = max I max diami^, max_ diami^ I , it follows that 

K 

i, (lj )= / fa (x) 2 dv g = Y. V 9 i F i)<t>i te) 2 + °(^) as v 0, 

/■ * 
7,K C ) = / ^(x) 2 dV 9 = £ V,(fi)0i(6) 2 + o(v d ) as 77 -j. 0, 

•^o i=l 



(60) 



for every j ^ jo- Note that, since is the complement of u>o in O, there holds 

Jj (w ) + ^ H) = / (z) 2 dV g + [ fa {xf dV g = l, (61) 



(62) 



for every j. Note also that 

k R 

Y,V g {Fi) = LV a {n) and = (1 - £)V^(fl). 

i=l i=l 

Set = (1 - 1)^(1*) and £ s = LV g {Fi). Then, we infer from (|60]) and (jUJ) that 

if 

(1 - I) 1,-M = Mi(6) 2 + o(v d ) as v -> 0, 

i=l 

II, (w ) =I-^^^-(6) 2 +o(?? d ) as 77 ->0, 
for every j jo- For e > to be chosen later, define the perturbation lj e of ujq by 

= U>\ (J ^(5^0 ) (J U 

V 1=1 / 8=1 

where = sh\ /n /Vg{B{&, l)) 1 /™ and £ = el\ ,n /V g {B{li, l)) 1 /™. Note that it is possible to define 
such a perturbation, provided that 

< £ < mm mm — — —, — — , mm - —. . 

h/ n l<i<* t\ jn ) 

It follows from the well known isodiametric ine quality and from a compactness argument that 
there exists a constant V n > (only depending on f2) such that V g (Fi) ^ ^(diam-Fj)™ for every i £ 
{!,••■ ,K}, and V g {Fj) ^ ^(diam i^)™ for every i E {!,••• , K}, independently on the partitions 



10 The isodiametric inequality states that, for every compact K of the Euclidean space R", there holds \K\ ^ 
|B(0,diam(A:)/2)|. 
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considered. Again, by compactness of fl, there exists v n > (only depending on fl) such that 
V g (B(x, 1)) ^ v n for every x £ CI. Set e = min(l, c$v n /Vn n )- Using (|5§)l . we get 

h\ /n " (l_L) 1 /-y n 1 /™diamF l " °' 

for every i G {1, — , K}, and similarly, 

7^(3(6, l)) 1 /" 



i/„ 



for every i G {1, •■• , It follows that the previous perturbation is well defined for every 
e £ (0,£o). Note that, by construction, 



K K 



V g {u s ) = V g {u a ) -^^(iJ&.l)) + Y d ?!V g {B{t h l)) 

i=l i=l 
K K 

= V g {u Q )-s n Y,hi + s n Y,Zi 



k k 



= V 9 {uj q ) £ ™(1 - L) £ V g (F t ) + E »L£ V fl (F,) 

i=l i=l 

= V 9 (lj ) - - L)LV g (fl) + e n L(l - L)V g (Cl) 

= v g (u ) = Lv g (n). 



Moreover, one has 



K „ K 



h = / h dV g = I s (wo) - X) / & (^) 2 dV g +J2 f 4>i dVJr 

and using again the fact that the ^ and are Lebesgue points of the functions (fyj, for every j ^ jo, 
we infer that 

k R 
i=i i=i 

= /i(wo) - e" f E Mi(6) 2 " E^te) 2 ) + °( J ? d ) as ^ ^ 0, 

and hence, using (1521) . 

1^) = Ijffjo) +e n (L- J,-(w )) + £ n o(r? d ) as 77 -> 0, 
for every j < j and every e g (0, £0). Since Eq ^ 1, it then follows that 

I 3 {u e ) ^ J(u ) + e n (L - J{uj q )) + e n o(r) d ) as rj -> 0, (63) 
for every j ^ jo and every e £ (0, £q), where the functional J is defined by pip . 
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We now choose the subdivisions (|58p fine enough (that is, rj > small enough) so that, for 
that 



every j ^ jo, the remainder term o (rj d ) in ([531 is bounded by h(L — J(wo)). It follows from (|63 



7^(o; E )> J( Wo ) + y(L-J( Wo )), (64) 
for every j ^ jo and every e <E (0, £o). 

Let us first show that the set ui e still satisfies an inequality of the type (|57l) for e small enough. 



Using (|54|) and Holder's inequality, we have 



\Ij(u e )- 7, Ml 



(Xw 0) - Xu (x)) <j>j (x) 2 dV g 

1/9 

<:A 2 ( / \x^(x)-x U0 (x)\ q dV g 
for every integer j and every £ € (0, £o), where g is defined by ^ + | = 1. Moreover, 

|;wO) -Xu (x)\ q dV g = / |x^(x) -XwoWl^ 

Jn 

(K K 
<=1 i=l 

= 2£ n L(l - L)V„(n), 
|/> £ ) - ^-(wo)! < (2A 2 V l L(l - L)F g (fi)) 1/9 

Therefore, setting 

(L - J(wo))« 



£i = mm ^£o, 
it follows from (|5"Tf that 



2 2 ?+i^4 2 ?i(l — L)V g {Vt) 



I J (Lo e )^L--{L~J{uj )) (65) 



1 

2 ( 

for every j ^ jo and every £ € (0, £i). 

Now, using the fact that J(uq) + %-(•£ — ^( w o)) ^ L — |(L — J(w )) for every £ € (0,£o), we 
infer from (|64| and (|65|) that 

£™ 

J(w e ) > J(«o) + —{L - JM), (66) 

for every e £ (0, £i). In particular, this inequality holds for e such that e n = C\ min(C2, L — J(o>o))> 
where the positive constants Ci and C2 are defined by 



8AL(1- iJV^fi) ' 2»Ci' 
For this specific value of £, we set wi = a; 6 , and hence we have obtained 

J(«i) ^ J(wo) + *y min(C 2 , L - J(w )) (£ - J(wo)). (67) 
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Note that the constants involved in this inequality depend only on L, A and CI. Note also that by 
construction wi satisfies a 5-cone property. 

If J(loi) ^ L then we are done. Otherwise, we apply all the previous arguments to this new 
set using QUE, there exists an integer still denoted jo such that (|57| holds with Uq replaced 
with u)±. This provides a lower bound for highfrequencies. The lower frequencies j ^ jo are then 
handled as previously, and we end up with (|64[) with lo q replaced with u)\. Finally, this leads to 
the existence of u>2 such that (|67[) holds with wi replaced with L02 and lu replaced with uj\ . 

By iteration, we construct a sequence of subsets uj^ of CI (satisfying a <5-cone property) of 
measure V g (ujk) — LV g (Cl), as long as J(wfc) < L, satisfying 

J{u)k+i) ^ J{uk) + — min(C 2 ,£ - J(wfe)) (L - J(w fe )). 

If J(uJk) < L for every integer /c, then clearly the sequence J(wfc) is increasing, bounded above by 
L, and converges to L. This finishes the proof. 

Remark 25. It can be noted that, in the above construction, the subsets are open, Lipschitz 
and of bounded perimeter. Hence, if the second problem is considered on the class of measurable 
subsets uj of CI, of measure V g (uj) = LV g (Cl), that are moreover either open with a Lipschitz 
boundary, or open with a bounded perimeter, then the conclusion holds as well that the supremum 
is equal to L. This proves the contents of Remark [T3l 



4.7 Proof of Proposition [2] 

Assume that CI is the unit (Euclidean) disk of 1R 2 , CI = {x £ 1R 2 | \\x\\ ^ 1}. It is well known that 
the normalized eigenfunctions of the Dirichlet-Laplacian are a triply indexed sequence given by 

rh ! m - / R oki r ) ir 3 = 0- 

pjhmW) - j Rjk{r)Yjm (e) if i ^ 1, 

for j G IN, k £ IN* and m — 1, 2, where (r, 0) are the usual polar coordinates. The functions Yj m (9) 
are defined by Yji(8) — ^= cos(j8) and Yj 2 {0) — sm(j6), and the functions iijfc are defined by 

where Jj is the Bessel function of the first kind of order j, and Zjk > is the fc th -zero of Jj. 
The eigenvalues of the Dirichlet-Laplacian are given by the double sequence of — z? k and are of 
multiplicity 1 if j = 0, and 2 if j ^ 1. Many properties are known on these functions and, in 
particular (see |40) ) : 

• for every j £ IN, the sequence of probability measures r 1— > Rjk(r) 2 rdr converges vaguely to 
1 as k tends to +00, 

• for every fc £ IN*, the sequence of probability measures r > Rjk(r) 2 rdr converges vaguely to 
the Dirac at r — 1 as j tends to +00. 

These convergence properties permit to identify certain quantum limits, the second property ac- 
counting for the well known phenomenon of whispering galleries. Less known is the convergence of 
the above sequence of measures when the ratio j/k is kept constant. Simple computations (due to 
[13j ) show that, when taking the limit of Rjk(r) 2 rdr with a fixed ratio j/k, and making this ratio 
vary, we obtain the family of probability measures 

1 r 

Ms = fs(r) dr = —=== X ( s ,i){r) dr, 
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parametrized by s £ [0, 1). We can even extend to s = 1 by denning /ii as the Dirac at r = 1. It 
easily follows that 



sup J(a) = sup inf / 
aeu L a eu L 3 m£ ' {1)2} J0 



2tt /.l 

2 



a{r,9)(j)jkm{r,9) rdrdd ^ sup K(a), 

aeu L 



where 

f 1 /-27T 



if(a) = inf / / \(r,9)d9f s (r) 



Lemma 7. There holds sup if (a) = £, and i/ie supremum is reached with the constant function 

a£U L 

a = L on Q. 

Proof of Lemma^ First, note that K(a = L) = L and that the infimum in the definition of K is 
then reached for every s € [0,1]. Since K is concave (as infimum of linear functions), it suffices 
to prove that (DK(a — L),h) ^ (directional derivative), for every function h defined on f2 such 
that Jjj h(x) dx = 0. Using Danskin's Theorem (see OE]), we have 

pi p2ir 

(DK(a = L),h)= inf / / h(r,9) d9 Ur) dr. 
se[o,i]J J 

By contradiction, let us assume that there exists a function h on f2 such that h{x) dx = and 
such that 

/i(r,0) d8f s {r) dr > 

/o Jo 

for every s £ [0,1]. Then, it follows that 

r-l r 2ir 



[ [ h{r, 9) d8 r dr>0 
Js Jo Vr 2 -s 2 



7T 

2 Jo 



Is JO 

for every s € [0, 1], and integrating in s over [0, 1], we get 

0< / / / h{r,9)d9^==drds= / —==ds h(r,9)d9dr 
Jo Js Jo Vr - s z Jo Jo Vr - s z Jo 

p2-K 

r I h(r,9)d9dr 
Jo 

7T f 

= — J h(x) dx = 0, 

which is a contradiction. The lemma is proved. □ 

It follows from this lemma that sup ag ^ J( a ) — L (note that a = L realizes the maximum), 
and hence, sup XujeUL J(Xu) ^ L. To prove the no-gap statement, we use particular (radial) subsets 
uj, of the form 

u = {(r,9) e [0,1] x [0,2tt] I 9euj e }, 

where \ujg\ = 2Lir, as drawn on Figure [5j 
For such a subset uj, one has 

<p jkm (xf dx = f R jk {rfrdr [ Y jm {6f d9 = f Y ]m (9) 2 d9, 



OJg J Log 
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A maximizing sequence for L-0.3 






Figure 5: Particular radial subsets 



for all j G N*, k G IN* and m = 1, 2. For j = 0, there holds 



(j>Qkm{x) 2 dx = I R jk {r) 2 rdr / d0 = |w e | 



Besides, since 



Ltt = |0| = / rdr / t»=i|w fl | > 

JO J mo A 



it follows that | org | = 2Ltt. By applying the no-gap result in dimension one (clearly, it can be 
applied as well with the cosine functions), one has 

sup inf / sin 2 (jf5) d8 = sup inf / cos 2 (jf5) d9 = Lit. 



\u> e \=2Ln 

Therefore, we deduce that 
and the conclusion follows. 



\oj e \=2Ln 



sup inf / (f>jk m {x) 2 dx = L, 



ie{i,2} 



4.8 An intrinsic spectral variant of the second problem 

The second problem (jTTJ) , defined in Section 11.11 depends a priori on the orthonormal Hilbcrtian 
basis (</>j)jeiN* of L 2 (fl) under consideration, at least whenever the spectrum of A is not simple. 
In this section we assume that the eigenvalues (A 2 )j 6 ]N* of A are multiple, so that the choice of 
the basis (<£j)ieJN* enters into play. 

We have already seen in Theorem [T] (see Section |2?2|) that, in the case of multiple eigenvalues, 
the spectral expression for the time-asymptotic observability constant is more intricate and it does 
not seem that our analysis can be adapted in an easy way to that case. 

Besides, recall that the criterion J defined by (jTTJ) has been motivated in Section 12721 by means 
of randomizing initial data, and has been interpreted as a randomized observability constant (see 
Thcorem[2]), but then this criterion depends a priori on the preliminary choice of the basis (^Jjgu* 
of cigenfunctions. 
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In order to get rid of this dependence, and to deal with a more intrinsic criterion, it makes sense 
to consider the infimum of the criteria J defined by over all possible choices of orthonormal 
bases of eigenfunctions. This leads us to consider the following intrinsic variant of our second 
problem. 

Intrinsic uniform optimal design problem. We investigate the problem of maxi- 
mizing the functional 

JintiXw) = inf / 4>{xfdV g , (68) 

over all possible subsets u> of SI of measure V g (oj) = LV g (Q), where £ denotes the set of 
all normalized eigenfunctions of A. 

Here, the word intrinsic means that this problem does not depend on the choice of the basis of 
eigenfunctions of A. 

As in Theorem [2j the quantity i-Jint(Xw) (resp., TJ; nt (xw)) can be interpreted as a con- 
stant for which the randomized observability inequality pip for the wave equation (resp., (|34l) 
for the Schrodinger equation) holds, but this constant is less than or equal to Axu) ( res P-i 

^randO^))- Besides, there obviously holds C^\xw) < |^int(Xw) and C^\xoj) < TJ int (Xu)- 
Indeed this inequality follows form the deterministic observability inequality applied to the partic- 
ular solution y(t, x) = e tXt (j)(x), for every eigenfunction <j). In brief, there holds 

Cf°(X") < l-M*-) < 4Tid(x-), and 4 5) ( X .) < rj int ( Xu ) < C ( T % ad ( Xu ). 
As in Section al! the convexified version of the above problem consists of maximizing the functional 

^int(a) = inf / a(x)(j)(xf dV„, 
over the set Ul- This problem obviously has at least one solution, and 



sup inf / \u{x)(l)(x) dV g ^ sup inf / a{x)4>{x) dV g . 

Theorem 6. Assume that the uniform measure via) ^9 * s a c ^ osure point of the family of prob- 
ability measures jj,^ — (p 2 dV g , <fi £ £ , for the vague topology, and that the whole family of eigen- 
functions in £ is uniformly bounded in L°°{VL). Then 

sup inf / (t>{xf dV g = sup inf / a(x)(j>(x) 2 dV g = L, (69) 

for every L 6 (0, 1). In other words, there is no gap between the intrinsic uniform optimal design 
problem and its convexified version. 



Proof. The proof follows the same lines as in Section 14751 by considering the function / defined by 
f{x) = (0(x) 2 ) 06£ . Then / <E L 1 ^, X) with X = L°°(£,R) which is a Banach manifold that can 
be seen as an infinite product of spheres of dimension equal to the respective multiplicities of the 
eigenvalues. □ 

Similarly, the intrinsic counterpart of Theorem [5] is the following. 
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Theorem 7. Assume that the uniform measure v dV g is the unique closure point of the family 

of probability measures /i^ = <fi 2 dV g , <fi G £ , for the vague topology, and that the whole family of 
eigenf unctions in E is uniformly bounded in L 2p (Q), for some p G (1, +00]. Then 

sup inf / <j){x) 2 dV g = L, (70) 

for every L G (0, 1). 

Proof. The proof follows the same lines as in Section 14.51 replacing the integer index j with the 
continuous index A (standing for the eigenvalues of ^4). The only thing that has to be noticed is 
the derivation of the estimate corresponding to (|64|) . In Section l4~5l to obtain (|64|) from ([63| . we 
used the fact that only a finite number of terms have to be considered. Now the number of terms 
is infinite, but however one has to consider all possible normalized eigcnfunctions associated with 
an eigenvalue |A| ^ |Ao|. Since this set is compact for every Ao, there is no difficulty to extend our 
previous proof. □ 

With respect to RemarkfTSl it is interesting to note that, here, we are able to provide examples 
where there is a gap between the intrinsic second problem (|68|) and its convexified version. 

Proposition 3. In any of the two following examples: 

• fl = S 2 , the unit sphere in M 3 , endowed with the usual flat metric; 

• f2 is the unit half-sphere in M 3 , endowed with the usual flat metric, and Dirichlet conditions 
are imposed on the great circle which is the boundary ofQ; 

if L is close enough to 1 then sup x ^ eWL J(Xu) < L, and hence there is a gap between the problem 
(|68[) and its convexified version. 

Proof. Assume first that f2 = S 2 , the unit sphere in H 3 . As mentioned in Remark [T9l in [37] it 
is proved that the set of semi-classical measures on S 2 coincides with the convex set of invariant 
probability measures for the geodesic flow that are time-reversal invariant. In particular, the Dirac 
measure /i 7 of any great circle 7 on S 2 (defined as an equator, up to a rotation) is the projection 
of a semi-classical measure. The measure /x 7 is the arc-length measure defined by 

If 1 

for every measurable subset u) of S 2 . Besides, since the uniform measure is a quantum limit as well, 
S 2 satisfies WQE and hence sup ng ^ J(a) = L (and the supremum is reached with the constant 
function a — L). Denoting by a the Lebesgue measure of S 2 , Ul is the set of all measurable subsets 
w of S 2 of measure a(L) = inL. For every oj G Ul, one has 

4ttL = / / Xw(<A 0) sm pdipdO ^ sine / / XwCv 3 ) 9) dipdO ^ sine / n uj\ — 2e)d9, 

J0 J0 Jo Je Jo 

for every e G [0, n/2], where 7^ denotes the great circle joining the north pole to the south pole at 
longitude 8 (where a north pole is fixed arbitrarily). By contradiction, assume that // 7e (w) > 3i/4 
for every G [0, 2tt]. Then we infer that A-kL > 27rsine(37rL/2 — 2e), which raises a contradiction 
when choosing e.g. e = n/ A and L close to 1. It then follows that 
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for every ui £ Ul, whence the gap. 

Assume now that Q is the unit half-sphere of H 3 . As recalled above, for every great circle 
C of S 2 there exists a sequence of squares of eigenfunctions fa whose support concentrates along 
C. Let S denote the orthogonal symmetry with respect to the hyperplane passing through the 
origin, cutting S 2 into two half-spheres, one of which being J7. Then, fa = (fa — fa o S)/^/2 is an 
eigenfunction of the Dirichlet-Laplacian on Q. Let us provj"] that the support of 'fa 2 concentrates 
on the union of two symmetric half-circles of fi, as drawn on Figure [51 



Figure 6: The half-sphere. 

Indeed, since fa 2 = \(4> 2 + \4>j S\ 2 + <j)j • fa °<S)> it suffices to prove that for every a € L°°(Q), 
J n a(x)4>j{x)4>j oS(x) da(x) tends to as j tends to +00. But this fact is obvious since the measure 
of the intersection of the corresponding supports tends to 0. The following (interesting in itself) 
fact follows: the Dirac measure along every union of symmetric half-circles on f2 is the projection of 
a semi-classical measure. Note however that, in this construction, the half-circles passing through 
the lowest point of the half-sphere cannot be considered. 

Then, the same calculation as before can be led. Indeed, let us fix a point N of the boundary of 
fl, and let S be the diametrically symmetric point, as on FigureEl If we think of N and S as a north 
pole and south pole, then any curve consisting of the union of two symmetric half-circles emerging 
from N and S can be viewed, with evident symmetries, as a great circle 7© of S 2 as considered 
previously. Then, the same argument can be applied and leads to the desired conclusion. □ 

5 Spectral approximation of the uniform optimal design 
problem 

In this section, we consider a spectral truncation of the functional J defined by (TTTj) , and we define 

Jn(Xu) = min fa (x) 2 dV g , (71) 

for every iV £ IN* and every measurable subset oj of f2, and we consider the spectral approximation 
of the second problem (uniform optimal design problem) 

sup J N {Xu>)- (72) 

11 This idea emerged from discussions with Luc Hillairet. 
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As before, the functional Jn is naturally extended to Ul by 

JnM = mm / a(x)6j(x) 2 dV a , 

for every a £ Ul- 

5.1 Existence, uniqueness and T-convergence properties 

One has the following result. 

Theorem 8. 1. For every measurable subset uj of CI, the sequence (Jn(Xu))ngin* is nonin- 
creasing and converges to J(Xuj)- 

2. There holds 

lim max Jn(o>) — max J (a). 

jV-H-oo aeUh aeUh 

Moreover, if (a N ) ne ]N* is a sequence of maximizers of Jn inU^, then up to a subsequence, 
it converges to a maximizer of J vtiUl for the weak star topology of L°° . 

3. Assume that M is an analytic Riemannian manifold and that Cl has a nontrivial boundary. 
Then, for every N £ IN*, the problem (1721) has a unique solution Xu N , where ui N £ Ul- 
Moreover, uj n is semi-analytic (see Footnote^) and has a finite number of connected com- 
ponents. 

Proof. For every measurable subset lu of Cl, the sequence (Jn (Xw))jveiN* is clearly nonincreasing 
and thus is convergent. Note that 




Hence, for every (ctj)j e w* £ £ 1 (M. + ), one has 

N . N 

1=1 Juj 1=1 

for every N £ IN* , and letting N tend to +oo yields 




and thus lim Jn(Xui) ^ J{Xu)- This proves the first item since there always holds Jn{Xu>) ^ 

TV— >+oo 

Since Jn is upper semi-continuous (and even continuous) for the L°° weak star topology and 
since Ul is compact for this topology, it follows that Jn has at least one maximizer a N £ Ul- Let 
a £Ul be a closure point of the sequence (a N ) ne w in the L°° weak star topology. One has, for 
every p $J N , 

sup J(a) ^ sup Jjv(a) = Jn(o- N ) ^ J p (a N ), 
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and letting N tend to +00 yields 



sup J (a) ^ lim Jjv(a ) ^ lim J P (a ) = J P {a), 

dfcM L 

for every p G IN*. Since J p (a) tends to J(a) ^ sup ag ^ L J(a>) as p tends to +00, it follows that a is 
a maximizer of J in W^. The second item is proved. 

To prove the third item, let us now prove that Jjv has a unique maximizer a" G Wi of Jn, 
which is moreover a characteristic function. We define the simplex set 

N 

An = {a= (aj)isy<iv I aj ^ 0, J^a, = 1}- 



Note that 



mm 



a(x)4>j (x) 2 dV g — min / a(x) ckj^ (x) 2 dV g , 



for every a <E Ul- It follows from Sion's minimax theorem (see [61]) that there exists a e .A 



A' 



such that (0^,0^) is a saddle point of the bilinear functional 

N 



(a,a)i-> / a(x) 7 a j<Pji x ) dV g 
Jn ~! 



2=1 

defined ouUl x -An, and 

iV iV 



max min / a(x) > aj(f>j(x) 2 dV g = min max / a(x) 7 aj(f)j(x) 2 dV c 
= max / a{x)Y]af ^{xf dV g = / a N {x) V af ^(x) 2 dV 9 . 



We claim that the function x i-> SjLi aI j 4>j{ x ) 2 IS never constant on any subset of positive 
measure. This fact is proved by contradiction. Indeed otherwise this function would be constant 
on VL (by analyticity) . At this step we have to distinguish between the different boundary conditions 
under consideration. For Neumann boundary conditions, we infer that A ff (5^- =1 a 1 ^ 4>j{x) 2 ) = 
on d£l (by continuity), and therefore X^jLi a f ' ^j4>j{ x ) 2 = on cT2, whence the contradiction since 
the coefficients are nonnegative for every j € {1, • • • , N} and Yjj=i a f = !• F° r the other 

boundary conditions, we infer that the function x 1— >• Yl^—i exf ' 4>j{ x ) 2 vanishes on f2, which is a 
contradiction. 

It follows from this fact and from (|75|) that there exists \ N > such that 

AT 



2=1 

otherwise, 



» = 

I 

for almost every x € fi. Hence there exists cj^ E Ul such that = Xu N ■ Since the eigenfunctions 
<j)j are analytic in fi (by analytic hypoellipticity) , it follows that w" is semi-analytic and has a 
finite number of connected components. □ 



4G 



Remark 26. Note that the third item of Theorem[5]can be seen as a generalization of [31] Theorem 
3.1] and [531 Theorem 3.1]. We have also provided a shorter proof. 

Remark 27. It is proved in [3TJ [54] that, in the one-dimensional case ft = [0,7r] with Dirichlet 
boundary conditions, the optimal set wjv maximizing Jjy is the union of N intervals concentrating 
around equidistant points and that ujn is actually the worst possible subset for the problem of 
maximizing Jjv+i- This is the spillover phenomenon. 



5.2 Numerical simulations: maximizing sequences 

We provide hereafter several numerical simulations based on the modal approximation described 
previously. 

Assume first that 17 = [0,7r] 2 , the Euclidean two-dimensional square. We consider Dirichlet or 
Neumann boundary conditions. In the Dirichlet case, the normalized eigenfimctions of A are 

2 

(f>j,k{xi,x 2 ) = -sin(jzi)sin(fcE 2 ), 

for every (xi,x 2 ) € [0,7r] 2 , and in the Neumann case the sine functions are replaced with cosine 
functions. Let N £ IN*. We use an interior point line search filter method to solve the spectral 
approximation of the second problem 

sup Jjv(Xw). 



where 



Jn(Xu>) = min / / Xu(xi,x 2 )4>j,k(xi,x 2 ) 2 dxidx 2 . 



Some results are provided on Figure [7] in the Dirichlet case, and on Figure [8] in the Neumann case. 

Assume now that fl = {x € H 2 | |x|| ^ 1}, the unit Euclidian disk of R 2 . We consider Dirichlet 
boundary conditions. The normalized eigenfimctions of the Dirichlet-Laplacian are a triply indexed 
sequence given by 

QjkmW) - | Rjk{r)Y]m{ e) if j > 1, 

for j € IN, k £ IN* and m = 1, 2, where (r, 9) are the usual polar coordinates. The functions Yj m (9) 
are defined by Yj\{6) — A= cos(jO) and Yj 2 (6) — sin(j^), and the functions Rjk are defined by 

where Jj is the Bessel function of the first kind of order j, and Zj% > is the fc th -zero of Jj. 
The eigenvalues of the Dirichlet-Laplacian are given by the double sequence of — z 2 fc and are 
of multiplicity 1 if j = 0, and 2 if j ^ 1. In Proposition [21 a no-gap result is stated in this 
case. Some simulations are provided on Figure We observe that optimal domains are radially 
symmetric. This is actually an immediate consequence of the uniqueness of a maximizer for 
the modal approximations problem stated in Theorem [8] and of the fact that f2 is itself radially 
symmetric. 
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Figure 7: On this figure, Q, = [0,7r] 2 . Line 1, from left to right: optimal domain (in green) in the 
Dirichlet case for N — 2 (4 eigenmodes) and L £ {0.2, 0.4, 0.6}. Line 2, from left to right: optimal 
domain (in green) for TV = 5 (25 eigenmodes) and L G {0.2,0.4,0.6}. Line 3, from left to right: 
optimal domain (in green) for N = 10 (100 eigenmodes) and L G {0.2,0.4,0.6}. Line 4, from left 
to right: optimal domain (in green) for N — 20 (400 eigenmodes) and L E {0.2, 0.4, 0.6} 



6 Further comments 

In this section, we first consider in Section [6TTI classes of subsets sharing compactness properties, 
in view of ensuring existence results for the second problem (|45|1 (uniform optimal design). In 
Section 16.21 we show how our results for the second problem can be extended to a natural variant 
of observability inequality for Neumann boundary conditions or in the boundaryless case. In 
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Figure 8: On this figure, Q = [0,7r] 2 . Line 1, from left to right: optimal domain (in green) in the 
Neumann case for N = 2 (4 eigenmodes) and L g {0.2, 0.4, 0.6}. Line 2, from left to right: optimal 
domain (in green) for N = 5 (25 eigenmodes) and L G {0.2,0.4,0.6}. Line 3, from left to right: 
optimal domain (in green) for N = 10 (100 eigenmodes) and L £ {0.2,0.4,0.6}. Line 4, from left 
to right: optimal domain (in green) for N — 20 (400 eigenmodes) and L E {0.2, 0.4, 0.6} 



Section [6.31 we study a variant of observability inequality for Dirichlet, mixed Dirichlet- Neumann 
and Robin boundary conditions, involving a H 1 norm, and we show that the criterion J of the 
second problem has to be slightly modified. We discover that the corresponding second problem has 
a unique solution (that is, an optimal set), in contrast with the previous results, whenever L is large 
enough. Finally, in Section lo\4"l we show how the problem of maximizing the observability constant 
is equivalent to an optimal design of a control problem and, namely, to that of controllability in 
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Figure 9: On this figure, fl = {x G H 2 | |x|| ^ 1} and L = 0.2. Line 1, from left to right: optimal 
domain (in green) in the Dirichlet case for N = 1 (1 eigenmode), N — 2 (4 eigenmodes) and 
N = 5 (25 eigenmodes). Line 2, from left to right: optimal domain (in green) for N = 10 (100 
eigenmodes) and N = 20 (400 eigenmodes) 



which solutions are driven to rest in final time by means of a suitable control function. 

6.1 Uniform optimal design on other classes of admissible domains 

According to Lemma [5J we know that, in the one-dimensional case, the second problem (1451) is 
ill-posed in the sense that it has no solution except for L = 1/2. In larger dimension, we expect 
a similar conclusion. One of the reasons is that the set Ul defined by (|44p is not compact for the 
usual topologies, as discussed in Remark [TUJ To overcome this difficulty, a possibility consists of 
defining a new class of admissible sets, Vl cWt, enjoying sufficient compactness properties and to 
replace the problem (1431) with 

sup J(xu)- (74) 

Of course, now, the extremal value is not necessarily the same since the class of admissible domains 
has been further restricted. 

To ensure the existence of a maximizer \u* °f Q74I) . it suffices to endow Vl with a topology, 
finer than the weak star topology of L°° , for which Vl is compact. Of course in this case, one has 

J(Xw) = ma x J(Xu>) < sup J(x«). 

This extra compactness property can be guaranteed by, for instance, considering some a > 0, and 
then any of the following possibles choices 

Vl = {Xu G Ul \ Pn(w) < a}, (75) 
where Pq(ui) is the relative perimeter of ui with respect to VL, 

Vl = {Xuj € U L \ ||xj|w(n) < «}, (76) 
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where || • is the BV(fl)-noim of all functions of bounded variations on Q (see for example 

0]), or 

Vl = {x^ € I w satisfies the 1/a-cone property}, (77) 

(see Section 14.51 footnote [5J . Naturally, the optimal set then depends on the bound a under 
consideration, and numerical simulations (not reported here) show that, as a tends to +00, the 
family of optimal sets behaves as the maximizing sequence built in Section [5] in particular the 
number of connected components grows as a is increasing. 

6.2 Further remarks for Neumann boundary conditions or in the bound- 
aryless case 

In the Neumann case, or in the case dfl = 0, there is a problem with the constants, as explained in 
Footnote [5J In this section, let us show that, if instead of considering the observability inequalities 
dU) and ([5]), we consider the inequalities 

4 W) (X")ll(yV)llw)^W) < / T / (\dty(t,x)\ 2 + \y(t,x)\ 2 ) dV g dt (78) 

JO Ju) 

in the case of the wave equation, and 

4 S) (x-)lly°ll 2 ff2 (o. C ) < f I (\dty(t,x)\ 2 + \y(t,x)\ 2 ) dV g dt (79) 

JO Jul 

in the case of the Schrodinger equation (see [621 Chapter 11] for a survey on these problems), 
then all results remain unchanged. Accordingly, the functional Gt formerly defined by (1101) is now 
replaced with 

G T (Xu,)= / / (\dty(t lX )\ 2 + \ y (t,x)\ 2 ) dV g dt. 

Jo Ju 

Indeed, consider initial data (y ,?/ 1 ) € 7? 1 (fi,C) x L 2 (£l, C). The corresponding solution y 
can still be expanded as (|12l) . except that now (4>j)j^* consists of the eigenfunctions of the 
Neumann-Laplacian or of the Laplacc-Bcltrami operator in the boundaryless case, associated with 
the eigenvalues (— A|)j 6 ]N*, with Ai = and cf>i which is constant, equal to 1/ ^V g (Vl). The relation 
(|14[) does not hold any more and is replaced with 

+00 

ll(yV)ll 2 H W) x W) = E ( 2A ?M 2 + 2A ^I 2 + K + b j\ 2 ) • ( 8 °) 

i=i 

Following Section l2~2"1 we define the time asymptotic observability constant C^\x^j) as the largest 
possible nonnegative constant for which the time asymptotic observability inequality 

c^Hx.My^y^fHW^LW) ^ £ Jj\d ty (t,x)\ 2 + \ y (t,x)\ 2 ) dv g dt (si) 

holds, for all (y° : y l ) G H (O, C) x L 2 (il,C). Similarly, we define the randomized observabil- 
ity constant Cj^l n a(Xu) as tne largest possible nonnegative constant for which the randomized 
observability inequality 

o^LMUv^y^fHHo^LW) ^ E [£ J O^MI 2 + M^)l 2 ) dv g dtj (82) 
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holds, for all y°(-) £ H 1 ^, C) and € L 2 (il, C), where y v is defined as before by 

The time asymptotic and randomized observability constants are defined accordingly for the 
Schrodinger equation. 

Theorem 9. Let u be a measurable subset o/fl. 

1. If the domain £1 is such that every eigenvalue of the Neumann- Laplacian is simple, then 

2C ) (x,) = ^(x.) = Ax*)- 



2. There holds 



Proof. Following the same lines as those in the proofs of Theorems Q] and [5] (see Sections 12.31 and 
we obtain 

4TL(xa,) = TcW( X a,) = Tr, 

with 

r = . nf Eg(l + A 2 )(K | 2 + l^f) I Mx? dV g 

((»i),(6,))^(C))"\{0} (2A?(|ojP + |6,| 2 ) + | aj + 6,| 2 ) 

Let us prove that T = ^J(xlu)- First of all, it is easy to see that, in the definition of T, it 
suffices to consider the infimum over real sequences (a,j) and (bj). Next, setting dj = Pj cos6j and 
bj = pjsmOj, since \aj + bj\ 2 = p 2 (l + sin(20j)), to reach the infimum one has to take dj = ir/4 
for every j g IN* . It finally follows that 



F = in i o£^ / <t>iW 2 dV g = -J( Xu ). 
( Pj )et 2 OR Ju 2 



+ 00 

1 



□ 



6.3 Variant of observability inequalities and optimal design results 

In this section we consider only the following boundary conditions: Dirichlet, mixed Dirichlet- 
Neumann (with To 7^ 0), Robin. We replace the observability inequalities ^ and ([5]) with the 
inequalities 

4 W) (x w )||(»°,tf 1 )||lrx(n,c)xi='{n,c)< / [ \d t y(t,x)\ 2 dV g dt (83) 

JO Jul 

in the case of the wave equation, and 

4 5) (x-)lly°ll^(ac) < f I \d t y{t lX )\ 2 dV g dt (84) 

JO Jui 

in the case of the Schrodinger equation. Note that here, even in the Dirichlet case, we consider 

the full H 1 norm defined by \\f\\m(n,c) = (ll/ll|a(n,c) + II ^ f\\h(n,c)) 1/2 ■ These inequalities 
hold true under GCC, as already mentioned. Of course these inequalities cannot hold in the 
boundaryless case or for Neumann boundary conditions, due to constants (as discussed previously). 
Since the norm used at the left-hand side is stronger, it follows that C^\xu>) ^ Cj^CXw) and 

4 s) (x*)<4 s W 
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Accordingly, the functional Gt formerly denned by (jlOp is now replaced with 

Gt(x^)= f I \d t y(t,x)\ 2 dV g dt. 

JO Juj 

In contrast to the previous results, the time averaging procedure or the randomization with 
respect to initial data does not lead to the functional J defined by (ITT]) but to the slightly different 
functional 

J(X.) = .inf. £ ^{x? dV g . 

We will see that the study of the associated optimization problem differs significantly from the one 
considered previously. 

Let us first explain how to derive this expression in the case of the wave equation (the arguments 
being similar for the Schrodinger equation). Consider initial data (y ,?/ 1 ) G i? 1 (fi, C) x L 2 (Q,C). 
One has 

+00 
i=i 

As previously, the time asymptotic observability constant C^\\ u ) is defined as the largest pos- 
sible nonnegative constant for which the time asymptotic observability inequality 



C t " ) (xMv°,V 1 )f„;a X )yLW> « T 1 ™ S / / \9,y(t,xfdV,dt (86) 



r 

holds, for all (y ,?/ 1 ) G ff 1 (fi,C) x L 2 (£1,C). Similarly, the randomized observability constant 
Crp^l nd (Xui) is the largest possible nonnegative constant for which the randomized observability 
inequality 

4Zd(X")ll(l/V)llWc)xLW) ^ E ^ Jjd tyi/ (t,x)\ 2 dV g dtj (87) 

holds, for all y°{-) G ^(0,0) and y x (-) G £ 2 (ft,C), where y„ is defined as before by (02]). The 
time asymptotic and randomized observability constants are defined accordingly for the Schrodinger 
equation. 

Theorem 10. Let lj be a measurable subset of ft. 

1. If the domain £1 is such that every eigenvalue of A is simple, then 

2CW(X.)=C< S) (X U ) = %). 

2. There holds 

Proof. Following the same lines as those in the proofs of Theorems [TJ and [2] (see Sections 12.31 and 
we obtain 



cgL(x«) = rcWOc w ) = rr, 

with 



T = inf 



((a 3 ),fe))e(^ 2 (C)) 2 \{o} (2A 2 (K-| 2 + |6,| 2 ) + | 0j + b 3 \ 2 ) 
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Let us prove that T = hJ(Xu)- First of all, it is easy to see that, in the definition of T, it 
suffices to consider the infimum over real sequences (a,j) and (bj). Next, setting a,j — pj cos9j and 
bj = pjSmOj, since \aj + bj\ 2 — p 2 (l + sin(20j-)), to reach the infimum one has to take dj = ir/4 
for every j g IN* . It finally follows that 



1 A z f 1 



□ 



We are thus led to introduce the following version of the second problem (uniform optimal 
design problem), associated with the above observability inequalities. 

Uniform optimal design problem. We investigate the problem of maximizing the 
functional 

J( Xu )= inf 7J / ^(xfdV g: (88) 

over all possible subsets uj of £1 of measure V g {ui) = LV g (Q), where the jj 's are defined 
by 

A? 

Note that the sequence ('Jj)jevi* is monotone increasing, and that < 71 ^ jj < 1 for every 
j G IN*. 

As in Section 14. 1[ the convexified version of this problem is defined accordingly by 

sup J(a), (89) 

a£U L 

where 

J(o)= .inf 7j / O(a0(fc(x) a d^. (9°) 
jew- J n 

As in Sections 14.11 and I4.2[ under the assumption that there exists a subsequence of (0^)j e ]N* 

converging to yjjfj m weak star L°° topology (i°°-WQE property), the problem has at least 

one solution, and sup ag ^- £ J{ a ) — L, and the supremum is reached with the constant function 
a = L. 

We will next establish a no-gap result, similar to Theorem |4l but only valuable for nonsmall 
values of L. Actually, we will show that the present situation differs significantly from the previous 
one, in the sense that, if 71 < L < 1 then the highfrequency modes do not play any role in the 
problem (|88]) . Before coming to that result, let us first define the truncated versions of the problem 
([55D . For every N € IN*, we define 

J N (a) = i m£ n l3 J a(x)^ (x) 2 dV g . (91) 

An immediate adaptation of the proof of Theorem [5] yields the following result. 
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Proposition 4. For every N £ IN* , the problem 

sup J N (a) (92) 
aeu L 

has a unique solution a N that is the characteristic function of a set lu n . Moreover, if M is analytic 
then uj n is semi-analytic and has a finite number of connected components. 

The main result of this section is the following. 

Theorem 11. Assume that the QUE and uniform L°° -boundedness properties hold. Let L £ 
(71, 1). Then there exists Nq £ IN* such that 

max J(xu) = max J N (Xu>) < 7i < L, (93) 



for every N ^ iVo- In particular, the second problem (|88[) has a unique solution Xui N o 7 an d moreover 
if M is analytic then the setuj No is semi-analytic and has a finite number of connected components. 

Proof. Using the same arguments as in Lemma 21 it is clear that the problem (|89|) has at least one 
solution, denoted by a°°. Let us first prove that there exists Nq £ INT* such that J(a°°) = Jn (o°°). 
Let e £ (0, L — 71). It follows from the L°°-QUE property that there exists Nq £ IN* such that 



7j / a°°(x)(j> j (xydVg>L-e, (94) 

for every j > Nq. Therefore, 



J(a°°) = j a°°(x)h( X ydV g 



^ min 



inf 7j / a°°{x)4> 3 {xYdV gi .inf 7j / a°°(x)4> 3 (x) 2 dV g 



(j No (a°°),L-e) = J No (a°°), 



since L — e > 71 and JN (a°°) ^71. It follows that J(a°°) = Jjv (a°°). 

Let us now prove that J(a°°) = J No (a N °), where o" is the unique maximizer of Jjv (see 
Proposition S]). By definition of a maximizer, one has J(a°°) — J No (a°°) < JN (a N °). By contra- 
diction, assume that JN {a°°) < JN (a N °). Let us then design an admissible perturbation at £Ul 
of a°° such that J(a*) > J(a°°), which raises a contradiction with the optimality of a 00 . For every 
£ S [0, 1], set a t = a 00 + t(a N ° — a°°). Since Jat is concave, one gets 

J No (a t ) t)J No (a°°) + tJ No (a N °) > J No (a°°), 

for every t £ (0, 1], which means that 

inf 7j / atix^ixf dV g > inf 7j / a°°{x)<t> 3 (xf dV g > J(a°°), (95) 

for every t £ (0, 1]. Besides, since a N ° — a°° £ (—2, 2) almost everywhere in tt, it follows from (|M|) 
that 

'' u N ° 

a 



7j / a t (x)(t> 2 j (x)dV g =-f j / a°°(x)^(x) 2 dy g + i 7j - / {a N ° (x) - a 00 (x))<j) 3 {x) 2 dV g ^ L - e - 2t, 
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for every j ^ N . Let us choose t such that < t < L | 71 , so that the previous inequality yields 

7j / ot(x)^(a;) 2 dK fl > 71 > 71 / a°° {x)<j>i{xf dV g > J(o°°). (96) 
Jn Jn 

for every j J? Ao- Combining the low modes estimate (|95p with the high modes estimate (I96[) . we 
conclude that 



J{a t ) = inf 7j / a^x)^)'^ > J(a°°), 

which contradicts the optimality of a°° . 

Therefore Jjv (a°°) = J(a°°) = J No (a Na ), and the result follows. □ 

Remark 28. Under the assumptions of the theorem, there is no gap between the second problem 
(|88p and its convexified formulation (|89|) . as well as before. But, contrarily to the previous results, 
here there always exists a maximizer in the class of characteristic functions whenever L is large 
enough, and moreover, this optimal set can be computed from a truncated formulation ([91]) for a 
certain value of N. In other words, the maximizing sequence {Xu N ) /vgn* resulting from Proposition 
S]is stationary. Here, the high modes play no role, whereas in the previous results all modes had the 
same impact. This is due to the fact that, in the left hand-side of the observability inequalities (1551) 
and (|84|) . we use the 7J 1 -norm instead of the norm of D(A 1 ^ 2 ). This causes a spectral asymmetry 
that finally leads to the above result. 

Remark 29. Here, if L is not too small then there exists an optimal set (sharing nice regularity 
properties) realizing the largest possible time asymptotic and randomized observability constants. 
The optimal value of these constants is known to be less than L but its exact value is not known. 
It is related to solving a finite dimensional numerical optimization problem. 

Remark 30. In the case where L ^ 71, we do not know whether or not there is a gap between 
the second problem (|B"8"]l and its convexified formulation (|8"9"]) . Adapting shrewdly the proof of 
Theorems H] or [S] does not seem to allow one to derive a no-gap result. Nevertheless, one can prove 
using these arguments that swp XuiGUL J(Xu) ^ JiL- 

Remark 31. We formulate the following two open questions. 

• Under the assumptions of Theorem [TT| does the conclusion hold true for every L 6 (0, 1)? 

• Does the statement of Theorem [TT1 still hold true under weaker ergodicity assumptions, for 
instance is it possible to weaken QUE into WQE? 

Remark 32. The QUE assumption made in Theorem [TTJ is very strong, as already discussed. 
It is true in the one-dimensional case but up to now no example of a multi-dimensional domain 
satisfying QUE is known. 

Anyway, we are able to prove that the conclusion of Theorem [TT] holds true in a domain which is 
a tensorized version of a one-dimensional domain. Indeed, consider either the domain fl = T" (flat 
torus), or the domain f2 = [0, 7r] n , the Euclidean n-dimensional square, with Dirichlet boundary 
conditions, or mixed Dirichlet-Neumann boundary conditions with either Dirichlet or Neumann 
condition on every full edge of the hypercube. The normalized eigenfunctions of A are then 

n 

<P]i... ln {x u ...,x n ) = W_<p 0k (x k ), 

k=l 

for all (ji, . . . ,j n ) € (IN*)", where the </>.,■ 's are the normalized eigenfunctions of the corresponding 
one-dimensional case, i.e., <f>j(x) — ^J^sin{i:x) or <fij(x) — \J~^ cos(7rx) for every x G [0,7r]. Ob- 
viously, fi does not satisfy QUE (nor QE), but satisfies WQE, and moreover the eigenfunctions 
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4>j 1 ...j n are uniformly bounded in L°°{Q). Let us prove however that the equality (|93l) holds. More 
precisely, one has the following result. 

Proposition 5. There exists L £ (0, 1) and N £ IN* such that 

max J(xu) = max Jjv(x w ), (97) 

for every L £ [L$, 1) and every N ^ iVo. 

Proof. The proof follows the same lines as the one of Theorem [TT] Nevertheless, the inequality 
(|94|) may not hold whenever QUE is not satisfied and has to be questioned. In the specific case 
under consideration, is replaced with the following assertion: there exists Nq £ IN* such that 
for e > 0, there exists Lo £ (0, 1) such that 

7ji-"in / a °°( x ) ( l ) ji...j n { x ) 2 dx ^ L — e, 

for every L £ [Lo, 1) and for all (ji, . . . ,j n ) £ (IN*)™ such that min(ji, . . . ,j n ) ^ A*o. This assertion 
indeed follows from the following general lemma. 

Lemma 8. Let p £ L°°(0, M+) be such that J Q p(x)dx > 0. Then 



P(x)<f>n...j n (x) 2 dx > (J p(x)dx 



inf / p(x)(j>i, i (xY dx ^ F^l ( / p(x) ) > 0, 

(i!,...,j„)e(w*)« 

where F(x) — —{x — sinx) /or every x £ [0, tt] and i 71 !™! = F o ■ ■ ■ o F (n times). 

This lemma itself easily follows from |55[ Lemma 6] (case n = 1) and from an induction 
argument. □ 

We end this section by providing several numerical simulations based on the modal approxi- 
mation of this problem for = [0,7r] 2 , the Euclidean square, with Dirichlet boundary conditions 
on dQ H ({x2 = 0} U {x2 = it} and Neumann boundary conditions on the rest of the boundary. 
Note that we are then in the framework of Remark [32j and hence the conclusion of Proposition [97] 
holds true. The normalized eigenfunctions of A are then 

2 

4>j,k(xi,x 2 ) = -sm(jx 1 )cos(kx 2 ), 

IT 

for all (xx, X2) £ [0, tt] 2 . The eigenvalue A^fe associated with the eigenfunction <f)j^ is Xj k — j 2 + k 2 . 
Let N £ IN* . As in Section 15.21 we use an interior point line search filter method to solve the 
spectral approximation of the second problem sup X[jgWi Jn(Xu)i with 



Jn(Xu) = min / / Xu>{xi, x 2 )4>j,k{x\, x 2 ) dx\dx 2 . 

l^j,k^Nj Q J 

Some numerical simulations are provided on Figures [TU] and QTJ On Figure [TU1 the optimal domains 
are represented for L £ {0.2, 0.4, 0.6}. In the three first cases, the number of connected components 
of the optimal set seems to increase with N. The numerical results provided in the case L — 0.9 
on Figure [TT] illustrate the conclusion of Proposition [5] showing clear evidence of the stationarity 
feature proved in this proposition. 
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Figure 10: On this figure, Q — [0, 7r] 2 , with mixed Dirichlet-Neumann boundary conditions. Line 1, 
from left to right: optimal domain (in green) for N = 2 (4 eigenmodes) and L £ {0.2, 0.4, 0.6}. Line 
2, from left to right: optimal domain (in green) for N = 5 (25 eigenmodes) and L £ {0.2, 0.4, 0.6}. 
Line 3, from left to right: optimal domain (in green) for N = 10 (100 eigenmodes) and L G 
{0.2,0.4,0.6}. Line 4, from left to right: optimal domain (in green) for N — 15 (225 eigenmodes) 
and L e {0.2, 0.4, 0.6} 



6.4 Optimal location of internal controllers for wave and Schrodinger 
equations 

In this section, we investigate the question of determining the shape and location of the control 
domain for wave or Schrodinger equations that minimizes the L 2 norm of the controllers realizing 
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Figure 11: On this figure, il — [0,7r] 2 , with mixed Dirichlet-Neumann boundary conditions, and 
L = 0.9. Line 1, from left to right: optimal domain (in green) for N £ {1, 2}. Line 2, from left to 
right: optimal domain (in green) for N G {5, 10} 



null controllability. In particular, we explain why this optimization problem is exactly equivalent 
to the problem of maximizing the observability constant. For the sake of simplicity, we will only 
deal with the wave equation, the Schrodinger case being easily adapted from that case. Also, 
without loss of generality we restrict ourselves to Dirichlet boundary conditions. 

Consider the internally controlled wave equation on fi with Dirichlet boundary conditions 

d tt y(t, x) - Ay(t, x) = h u (t, x), (t, x) G (0, T) x O, 

y(t, x) = 0, (t, x) G [0, T] x 90, (98) 

2/(0, x) = y°(x), d t y(0, x) ^y 1 {x) : x£ fi, 

where h u is a control supported in [0,T] x uj and w is a measurable subset of il. Note that the 
Cauchy problem ([9"5]) is well posed for all initial data (y ,!/ 1 ) £ Hq(D,,C) x L 2 (il,C) and every 
K £ L 2 ((0,T) x Q,C), and its solution y belongs to C°(0, T; Hg(ft, C)) n C 1 (0, T; L 2 (f2, C)) n 
C 2 (0, T; i/ _1 (i7, C)). The exact null controllability problem settled in these spaces consists of 
finding a control steering the control system (l9"5)l to 

y(T,-) = d t y(T, ■) = (). (99) 

It is well known that, for every subset w of fl of positive measure, the exact null controllability 
problem is by duality equivalent to the fact that the observability inequality 

C\\(4> ,^)fL W)xH - W )< [ I \ct>(t,x)\ 2 dV g dt, (100) 

JO J u) 

holds, for all (</>°, cf) 1 ) £ L 2 (£l, C) x if _1 (0, C), for a positive constant C (only depending on T and 
w), where tf> is the (unique) solution of the adjoint system 

d tt (f>(t, x) - Acj)(t, x) = 0, (t, x) £ (0, T) x n, 

<f>(t, x) = 0, (t, x) £ [0, T] x dn, (101) 

<p(0,x)=4>°(x), dt</)(0,x) =<i> x (x), x£fl. 
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The Hilbert Uniqueness Method (HUM, see [UJI3B]) provides a way to design the unique control 
solving the control problem and having moreover a minimal L 2 ((0,T) x f2,C) norm. 

This control is referred to as the HUM control and is characterized as follows. Define the HUM 
functional J u by 

Ju ^ ^ = l£ J ^x) 2 dV g dt~{4,\y°) H - 1Ml +{4P,y l ) L .. (102) 

The notation (■, ^h-x^h 1 stands for the duality bracket between if _1 (fi, C) and Hq (f2, C), and the 
notation (•, -)l 2 stands for the usual scalar product of L 2 (Cl,C). If (|100[) holds then the functional 
J u has a unique minimizer (still denoted ((p ,^ 1 )) in the space £ 2 (f2, C) x i7 _1 (57, C), for all 
(y ,?/ 1 ) G H^(Q.,C) x L 2 (Q,,C)- The HUM control h u steering (j/ ,^ 1 ) to (0,0) in time T is then 
given by 

hu(t,x) = Xu(x)4>(t,x), (103) 

for almost all (t, x) G (0, T) x f2, where 4> is the solution of (|101[) with initial data (</>°, t/) 1 ) minimizing 
Jul- 

The HUM operator r w is defined by 

: Hi(n,c) * L 2 (n,c) — -> £ 2 ((o,t) x n,c) 

Optimal design control problem. We investigate the problem of minimizing the 
norm of the operator 

llr.H = sup { ll^ll^((o.T)xn, C ) 0; x E 1(fi; x m 1 

( 11(2/ >rJ \\H£(n,c) xL2(o,c) J 

(104) 

over the set Ul ■ 

Here, we formulate the optimal design control problem in terms of minimization of the operator 
norm of in order to discard the dependence with respect to the initial data (y , y 1 ) and improve 
the robustness of the cost function. The next result establishes that the problems (fTTj) and (|104p 
are equivalent, and hence that the approach developed in Sections [2] and [4] is also well adapted to 
this optimal design control problem. Simultaneously, we generalize |55j where similar issues were 
investigated in the one-dimensional case. 

Proposition 6. Let T > and let lu be measurable subset of CI. If Cy (Xw) > then 
and if C^p W \xuj) — 0, then \\T U \\ = +oo. 

Proof. Denote by <j) u the adjoint state solution of (|101[) whose initial data minimize the functional 
J u . Then <f>^ can be expanded as 

i=i 

where the sequences A = (A^)j S ]N* and B = {B")j^» belong to £ 2 (C) and are determined in 
function of the initial data (0°,</>i) minimizing J w . Since J w is convex, the first-order optimality 
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conditions for the problem of minimizing J w over L 2 (O, C) x H 1 (f2, C) are necessary and sufficient. 
In terms of Fourier coefficients, they are written as 

A U {A,B)=C, (105) 

where the operator A w : (£ 2 (C)) 2 -> (£ 2 (C)) 2 is defined by 

MA-B)i = ^ J^A^ + B k e-^)<p k (x)<P 3 {x)^ *_ iXjt jdV g dt, 

for every j € IN*, with the notation A U1 (A,B) = (A UJ (A, B)j)jew* , and where 



-(^•^ 1 )l 2 ,i, 2 



for every j € IN*. For all (A, B) e [^ 2 (C)] 2 , one has 



(A B (^B) ) (A,B)) ( p (C ) )J = 
and it follows that 



T 

Ju> 



fe=i 



^W, , . (A,(i,5),(A,B))( P( c)) 2 
T U«J ^ 11/4 5\p ^ ^ J • 

||(A i? )ll(^ (c)) 2 

Indeed, we obtain the left-hand side inequality by definition of the observability constant. The 
right-hand side one is easily obtained, writing that the integral of a nonnegative function over 
w is lower than the integral of the same function over f2, which permits to use the orthogonality 
properties of the </>j's. By duality, we deduce that A w is a continuous symmetric invertible operator 
from (£ 2 (C)) 2 to (£ 2 (C)) 2 . Note that 

Mr , n ( A J 1 (C); C '>(£ 2 (]R)) 2 II A " / ( C ')ll(i? 2 (R)) 2 

ll r -ll = SU P W7^2 = SUp — p , 

Ce(£ 2 (K.)) 2 \{0} ll°ll(£ 2 (]R)) 2 Ce(£ 2 (R)) 2 \{0} ll C ll(£ 2 (]R)) 2 

where A^ 1 ^ 2 denotes the square root of the operator A" 1 . Setting ip = A W 1//2 (C), one computes 

2 1 

l(<? 2 (R)) 2 1 

sup 



^ 2 (R)) 2 \ { o } ||Ai%)|| 2 , 2(R))2 . nf f^J^w 1 ^ e ( p(R )n{0} } 

1 1 



inf( <A ,,y ( - 2W)2 I ^ (^ 2 (IR)) 2 \{0}1 4 W) (X.)' 



IIVll (( ,2 ( i R)) 2 



The conclusion follows. □ 
It follows from this result that, for the optimal design control problem, 



inf ||rj= ( sup Cf°(xj) 1 



and therefore the problem is equivalent to the problem of maximizing the observability constant. 
Then, all considerations done in this article can be applied to the optimal design control problem 
as well. 
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A Proof of Proposition [T] 



We focus first on the one-dimensional wave equation. Whereas the proposition is stated on [0, 7r], 
we assume hereafter that we are on [—it, it] , in order to facilitate the use of Fourier series. To avoid 
any technical problem we assume that we are in the framework of Remark [3] In particular, we 
assume that T is an integer multiple of 27r. According to the characterization of the optimal set 
in terms of a level set of the function 

+00 

3=1 

with atjj — pTr(aj + tij), and noting that the coefficients otjj are nonnegative and of converging 
sum, it suffices to prove the following result. 

Proposition 7. There exist a measurable open subset C of [— 7r,7r], of Lebesgue measure \C\ £ 
(0, 2tt), and a smooth function f on [— n,ir], satisfying the following properties: 

• C is of fractal type, and in particular has an infinite number of connected components; 

• f{x) > for every x £ C, and f(x) = for every x £ [— 7r, 7r] \ C; 

• f is even; 

• for every integer n, 

a n = / f(x) cos(ncc) dx > 0; 

J — 7T 

• The series ^2 a n is convergent. 

Proposition [TJ follows from that result, and the optimal set u> is then the complement of the 
fractal set C. By considering cartesian products of this one-dimensional fractal set, it is immediate 
to generalize the construction to a n-dimensional hypercube for the Schrodinger equation since the 
solution remains periodic in this case, which ensures that Gt does not involve any crossed terms. 

There are many possible variants of such a construction. We provide hereafter one possible 
way of proving this result. 

Proof. Let a € (0, 1/3). We assume that a is a rational number, that is, a — ^ where p and q 
are relatively prime integers, and moreover we assume that p + q is even. Let us first construct 
the fractal set C c [— 7T, it]. Since C will be symmetric with respect to 0, we describe below the 
construction of On [0, %]. Set s — and 

Sk =tt- + , 



for every k £ IN*. Around every such point Sk, k € IN*, we define the interval 

h 



of length |4| = 5 ^ T a(l-«) fc . 

Lemma 9. We have the following properties: 

• inf 1 1 > air; 

• sup/fc < inf 4+i < TT for every k £ IN* . 
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Proof. Since a < 1/3 it follows that inf I\ = tt — j(a + 1) > an. For the second property, note 
that the inequality sup 1^ < inf Ik+i is equivalent to 

a(l-a)* -1 (3-or) < (a + l) fc , 

which holds true for every k G N* since a(3 — a) < a + 1. □ 

It follows in particular from that lemma that the intervals Ik are two by two disjoint. Now, we 
define the set C by 

Cn[0,7r] = [0,a7r]U |J I h . 

fc=i 

The resulting set C (symmetric with respect to 0) is then of fractal type and has an infinite number 
of connected components (see Figure [T^)) . 

We now define the function / such that / is continuous, piecewise affine, equal to outside C, 
and such that f(sk) = bk for every fc 6 IN, where the bk are positive real numbers to be chosen 
(see Figure \T^i . 

Let us compute the Fourier series of /. Since / is even, its sine coefficients are all equal to 0. 
In order to compute its cosine coefficients, we will use the following result. 

Lemma 10. Let a G M, £ > and b > 0. Let g be the function defined on R by 

!^r(x — a+f) if a— f < £ < a, 
otherwise. 
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In other words, 



s, g is a positive triangle of height b above the interval [a — §]• Then 

f , . , , , 46 , . / raT 
/ g(x) coslna;) aa; = cos(na) 1 — cos — 
Jm ^ 2 \ 2 J 

for every n € IN* . 

It follows from this lemma that 



fix) cos(nx) dx = — — r(l — cosfncra")) 



airn z 



and 



f 2 k+1 bk ( tvk , / (nix . ■ 

f{x)cos{nx)dx = — _ a ^ fc7rn2 cos ^n7r - ^(a + 1) J (1 - cos - a) 



/ j a(l — a) k im 2 

for every fc e IN*. Note that 



/(x) cos(ncc) (ice 



46, 



avm 2 V 1 



(106) 



(107) 



(108) 



for 



every k £ IN* . Formally, the n th cosine Fourier coefficient of / is given by 

/7r pair +00 _ 

f (x) cos(nx) dx = 2 / /(a;) cos(na;) da; + 2 / f (x) cos(nx) dx . 
-7T JO fc=l"' / '= 



Our next task consists of choosing adequately the positive real numbers 6^, k 6 IN, so that the 
series appearing in the above formal expression of a n 
series of eeneral term a~ is convergent. 



iv^ciji iiuiiiiuvji ij *Jfc , nj ^_ Ji. 1 , until iiiv 

is convergent, a„ is nonnegative, and the 

series of general term a n is convergent. 

Let us first consider the integral (|106|) (first peak). It is clearly nonnegative for 
and is positive except whenever ??. is a. multinle r r 11 > i- ' • <■ •> • 



Let us hrst consider tile integral l|lUt)|) (hrst peak j . It is clearly nonnegative tor every n <E IN* , 
and is positive except whenever n is a multiple of 2q. Taking advantage of the rationality of a, we 
can moreover derive an estimate from below, as follows. Set 

(To = rninjl — cos(ti-vt) | n = 1, . . . , 2q — 1}. 
One has (Jq > 0, and there holds 



f(x) cos(nx) dx - ° , 
cmn 



for every n € IN* \ (2g]N*). At this step, assume that 



b k € 



1- a \ h 1 o-pbo 



2 ) 2 k 8 ' 



(109) 



(110) 



for every k £ IN* (60 > is arbitrary). Under this assumption, using (|108p it follows that the 
formal expression of a n above is well defined, and that 



+00 . 

/ f(x) cos (nx) dx 

k=l 



^ 2 ^ o / f(x)cos(nx)dx, 

2 airn z 2 J 



for 



every n £ IN* \ (2glN*), ensuring therefore a n > for such integers n. 
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If n — 2rq, with r £ IN*, then the integral (|106|) vanishes. We then focus on the second peak, 
that is, on the integral (|107[) with k = 1. Since n = 2rq, its value is 

/ fix) cosinx) dx = —, — ~~~r 7t cos ( 2rqn — rqiri— + 1) ) [ 1 — cos f rqw-tl — -) ] ] . 

7 7l a(l-a)7rn 2 \ q J \ \ q q J J 

Since p + g is even, it follows that cos (2rqn — rq-ni^ + 1)^ = 1. Hence, we have 

/ f(x) cosinx) dx = — — — ^ tt f 1 — cos ( rn-(q — pj \ ) ^ 0. 

J h a(l - a)im 2 \ \ q J J 

Moreover, since the integers p and q are relatively prime integers and q — p is even, in this last 
expression one has cos(m^(q — p)) = 1 if and only if r is multiple of q, that is, if and only if n is 
multiple of 2q 2 . As before we derive an estimate from below, setting 

<7i = min 1 1 — cos ^rir— (q — p)^j | r = 1, . . . , 2q — 1 

One has o\ > 0, and there holds 

/ f(x)cos(nx)dx ^ ( ^ hia \ — o' ( m ) 

for every n £ (2gN*) \ (2q 2 K*). At this step, additionally to pTU|) assume that 

1 — a \ 1 



for every fc ^ 2. Under this assumption, using (|108[) it follows that 



/ f( x ) cos(nx) dx 

l — o •> Ik 



^ \ i^ bl(J \ — 2 ^ \ [ f(x)cos(nx)dx, 
2 a(l — a)im z 2 J Il 



for every n 6 (2<7JM*) \ (2g 2 lN*), ensuring therefore a n > for such integers n. 

The construction can be easily iterated. At iteration m, assume that n = 2rq m , with r g IN*. 
Then the integrals over the m first peaks vanish, that is, 

f{x) cosinx) dx = f(x) cos(nx) dx = 

for every k = 1, . . . , m — 1. We then focus on the (m + l) th peak, that is, on the integral (jl07[) 
with k ~ m. Since n — 2rq m , its value is 

fix) cosinx ) dx = — r 7 cos 2rn 7r r — hi 

JK ' K ' a{\~a) m nn 2 \ H 2 m ^ \q 



rq-irp I p 
XI — cos r - 1 



Since p + q is even, it follows that 



•"^nfV-^fj + l) 1 = 1. 
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and hence, 



f{x) cos{nx) dx 



2 m+1 h 



a(l — a) m 7rn 2 



T7T p 



(<? - pY 



> o. 



Moreover, since the integers p and q are relatively prime integers and q — p is even, it follows easily 
that q and ( 2 ^ £ ) m are relatively prime integers, and therefore this last expression vanishes if and 
only if r is multiple of q, that is, if and only if n is multiple of 2q m+1 . Setting 

{/ rir p 

f(x) cos(nx) dx 

Ji m 

fi+- 

assumptions, we assume that 



one has a m > and 



r = l,...,2<z-lj>, 
2 m+1 b a 



J i a (l — a) m im 2 ' 

for every n £ (2q m fi*) \ (2q rn+1 ¥S*). Additionally to (|TTU1) . (ITT^|) and the following iterative 



b k ^ 



1 - a 



k-r 



1 



2fc— m+2 



(113) 



for every k ^ m + 1. Under this assumption, using (|108[) it follows that 



/J / f( x ) cos(nx) dV g 

fc=m+l 



< 1 2 m+1 b m a m < 1 
^ 2 a(l - a)" 1 ™ 2 ^ 2 



/(x) cos(nx) dV 9 , 



for every n 6 (2g" l ]N*) \ (2g" l+1 ]N*), ensuring therefore a„ > for such integers n. 

The construction of the function / goes in such a way by iteration. By construction, its Fourier 
cosine coefficients a n are positive, and moreover, the series J2n=o a « i s convergent. We have thus 
constructed a function / satisfying all requirements of the statement except the fact that / is 
smooth. 

Let us finally show that, using appropriate convolutions, we can modify / in order to obtain a 
smooth function keeping all required properties. Set fa = f[- a ir,a-n-] an d fk = fi k for every k £ INT*. 
For every e > 0, let p e be a real nonnegative function which is even, whose support is [— e,e], 
whose integral over H is equal to 1, and whose Fourier (cosine) coefficients are all positive. Such 
a function clearly exists. Indeed, only the last property is not usual, but to ensure this Fourier 
property it suffices to consider the convolution of any usual bump function with itself. Then, for 
every k G IN, consider the (nonnegative) function fk defined by the convolution fk = p E (k) * fk, 
where each e{k) is chosen small enough so that the supports of all functions fk are still disjoint 
two by two and contained in [—it, tt] as in Lemma |H1 Then, we define the function / as the sum of 
all functions and we symmetrize it with respect to 0. Clearly, every Fourier (cosine) coefficient 
of / is the sum of the Fourier (cosine) coefficients of fk, and thus is positive, and their sum is 
still convergent. The function / is smooth and satisfies all requirements of the statement of the 
proposition. This ends the proof. □ 
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